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ABSTRACT 


An  approximate  mathematical  model  was  formulated  to 
describe  the  "push-pull"  method  of  steam  injection  during 
the  backflow  period.,  Two  example  calculations,  in  which  the 
method  of  Willman  et  al(23)  was  used  to  calculate  the  in¬ 
jection  period, are  included. 

The  thermal  behavior  of  the  strata  surrounding  the 
permeable  sand  was  approximated  by  the  one-dimensional  heat 
flow  equation.  The  overburden  and  underburden  were  con¬ 
sidered  as  two  sets  of  concentric,  semi-infinite  cylinders, 
and  the  permeable  sand  was  assumed  to  have  an  infinite  ther¬ 
mal  conductivity  in  the  vertical  direction. 

The  temperature  distribution  in  the  sand  was  deter¬ 
mined  by  numerically  solving  the  total  energy  balance 
equation  for  the  system.  Darcy's  law,  employing  plug  flow, 
was  used  to  compute  the  saturation  distribution  and  composi¬ 
tion  of  the  fluid  flowing  at  various  time  steps  of  the 
numerical  solution. 

With  the  exception  of  viscosity,  all  physical  proper¬ 
ties  were  assumed  to  be  independent  of  temperature  and  were 
calculated  at  an  average  temperature. 

Although  a  verification  of  the  model  was  not  possible, 
the  general  shape  of  the  calculated  back  production  versus 
time  curves  is  similar  to  those  observed  in  field  operations (2) . 
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INTRODUCTION 

Statistics,  reported  by  Boberg(2),  indicate  that,  at 
the  end  of  1964,  60%  of  the  active  new-field  projects  for 
secondary  recovery,  in  the  United  States,  involved  thermal  pro¬ 
cesses,  as  compared  with  about  20%  in  1958.  This  increased 
popularity  of  thermal  recovery  processes  is  primarily  due  to 
higher  recovery  of  high  viscosity  oils  obtained  by  reservoir 
heating. 

The  various  thermal  techniques  proposed  may  be  divided 
into  four  basic  groups. 

1.  Electrical  methods,  in  which  an  electrical  heater  is 
installed  at  the  bottom  of  the  hole  to  heat  either  the 
fluids  being  produced  or  the  fluids  being  injected. 

2.  In  situ  combustion  projects,  in  which  part  of  the  oil-in¬ 
place  is  ignited  to  heat  the  reservoir. 

3.  Hot-fluid  injection  methods,  in  which  steam,  hot  water, 
hot  air  or  hot  gases  are  injected  into  the  reservoir. 

4.  Nuclear  explosion  techniques. 

The  injection  of  hot  fluids  combines  both  flooding  and 
thermal  effects.  These  thermal  effects  have  been  identified 
as  thermal  expansion  of  the  oil  and  viscosity  reduction (23) . 
Steam  supplies  a  large  heat  of  vaporization  and,  in  addition, 
may  cause  steam  distillation  which  further  increases  the  re¬ 
covery  of  oils  with  large  fractions  of  steam-distillable  com¬ 
ponents  ( 23)  . 
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The  major  disadvantage  of  thermal  projects  is  the 
large  heat  loss  to  the  strata  surrounding  the  permeable  reser¬ 
voir.  For  example,  calculations  of  the  thermal  efficiency, 
using  the  method  of  Marx  et  al(16),  indicate  that  the  heat 
losses,  after  five  months  of  injection  into  a  10-feet  thick, 
circular  reservoir,  may  be  as  high  as  80%  of  the  heat  injected. 
Further  analysis  shows  that  most  of  the  heat  present  in  the 
steam  and  hot  water  portions  of  the  sand  does  not  reach  the 
oil  bank  ahead  of  the  intermediate  cold  water  zone.  This  be¬ 
havior  of  the  injection  process  led  to  the  idea  of  reversing 
the  direction  of  flow  by  stopping  injection  and  producing 
through  the  injection  well,  to  recover  some  of  the  heat  stored 
in  the  sand.  In  addition,  if  conditions  are  such  that  the 
rate  of  back  production  is  high,  the  reservoir  gand  may  be¬ 
have  as  a  heat  sink  with  respect  to  the  surrounding  strata. 

This  alternating  process  of  injection  of  steam  and  production 
of  oil  has  been  called  the  "push-pull"  or  "huf f-and-puf f " 
method  of  steam  injection. 

To  date,  no  theoretical  study  of  this  process  has 
been  published.  The  purpose  of  this  work  is  to  present  «an 

i 

approximate  approach  to  the  problem  of  predicting  the  back 


flow. 
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THEORY  AND  LITERATURE  SURVEY 

Heat  Losses 

Saturated  steam  is  injected  into  a  horizontal,  homo¬ 
geneous,  isotropic  oil  reservoir  of  constant  thickness  (see 
Figure  la).  The  heating  of  the  flooded  sand  and  the  loss  of 
heat  to  the  surrounding  strata  causes  steam  condensation  and 
gradual  cooling  of  the  hot  water.  Three  distinct  thermal  zones 
may  be  visualized:  a  water  zone  at  reservoir  temperature,  a 
hot  water  zone  at  temperatures  varying  from  reservoir  to  steam 

temperature  and  the  steam  zone  at  the  saturated  steam  temperature. 

* 

The  temperature  distribution  in  the  flooded  sand  has 
the  shape  represented  by  the  continuous  line  in  Figure  1C. 

This  representation  neglects  the  variation  of  steam  temperature 
with  distance  in  the  steam  zone,  due  to  a  decrease  in  pressure 
with  an  increase  in  distance  from  the  point  of  injection. 

To  calculate  the  position  of  the  various  zones  at  any 

p- 

given  time,  it  is  necessary  to  take  into  account  the  following 
effects : 

1.  The  wellbore  heat  losses. 

2.  The  heat  losses  to  overburden  and  underburden. 

3.  The  heat  required  to  heat  the  flooded  sand. 

An  overall  heat  balance  for  the  system  may  be  written 
as 

Rate  of  heat  Rate  of  heat  Rate  of 

injection  at  =  injection  at  -  wellbore 

sand  face  surface  heat  losses 


Heating  rate 
+  of  flooded 
sand 


Rate  of  heat  losses 
to  overburden  and 
underburden 
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DIAGRAM  OF  THE  PHYSICAL  SYSTEM 


FIG.  I 
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Hf(t)  =  H.(t)  -  <Hl)wb 

=  <V0B  +  H(t)  U) 

Marx  et  al(16)  presented  an  analytical  solution  to  the 

problem,  employing  the  following  assumptions: 

lo  The  heat  injection  rate  at  the  sand  face  is  maintained 
constant o 

2.  A  negligible  temperature  field  is  created  in  the  overburden 
by  the  wellbore  heat  losses.  Thus,  the  underburden  has  the 
same  thermal  behavior  of  the  overburden. 

9 

3. j£  The  sand  and  the  surrounding  impermeable  strata  have  zero 

thermal  conductivity  in  the  horizontal  direction. 

4.  The  permeable  sand  has  an  infinite  thermal  conductivity  in 
the  vertical  direction. 

5.  The  actual  temperature  profile  in  the  permeable  sand  may  be 
approximated  by  the  dashed  line  shown  in  Figure  1C. 

6.  Radial,  one-dimensional,  incompressible  fluid  flow  is  tak¬ 
ing  place  in  the  reservoir. 

The  solutions  obtained  are 
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(Hl)ob  =  Hf  (1  -  e’jZ  erfc  u)  (5) 

The  radii  for  the  steam  and  hot  water  zones  may  be  calculated 
from  the  "effective  radius" 

Equation  2  was  modified  by  Willman  et  al(23)  to  estab¬ 
lish  formulations  for  the  saturation  distribution  in  the  sand 
during  the  steam  injection  process. 

Lauwerier (11)  obtained  an  analytical  solution  for  the 
hot  water  injection  case  with  the  same  assumptions  made  by 
Marx  et  al.  Spillette (22)  examined  contributions ,  in  Russian 
literature,  of  Avdonin(l)  and  Rubinshtein { 20} ,  for  hot  water 
injection,  and  solved  a  general  model  employing  assumptions  1 
and  6  only.  The  results  of  these  studies,  dealing  with  hot 
water  injection,  are  applicable  to  the  analysis  of  the  hot 
water  zone  advancing  ahead  of  the  steam. 

The  Controlling  Heat  Transfer  Mechanism  in  Porous  Media  with 

A  Flowing  Fluid 

The  theoretical  prediction  of  the  temperature  distri¬ 
bution  in  the  permeable  reservoir  sand  during  injection  or  back 
production  is  simplified  if  controlling  heat  transfer  mechanisms 
can  be  selected.  Two  physical  models  were  studied  for  the 
selection  of  the  predominant  mechanisms:  1)  a  semi-infinite 
porous  body  through  which  hot  water  flows  in  a  linear  direction, 
and,  2)  a  cylindrical  model  in  which  the  hot  water  injected  at 
r  =  0  flows  in  the  radial  direction.  These  two  models  are  re¬ 
presented  in  Figure  2.  The  solid  and  fluid  temperatures  in 
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both  models  are  initially  equal  to  the  same  constant  value 
throughout.  At  time  zero,  the  injected  water  temperature  is 
suddenly  changed  to  a  higher  value  and  is  held  constant  there¬ 
after.  The  problem  is  to  describe  the  fluid  and  solid  tem¬ 
peratures  as  a  function  of  time  and  position  from  the  point 
of  injection. 

In  a  porous  medium,  through  which  fluid  is  passing, 
thermal  energy  can  be  transferred  from  one  part  of  the  system 
to  another  by  four  basic  mechanisms (5) : 

1.  the  physical  movement  of  the  fluid  which  carries  its  own 
thermal  energy; 

2.  the  conduction  of  heat  through  the  solid  and  fluid  phases? 

3.  the  convective  heat  transfer  between  the  solid  and  fluid 
phases ,  and 

4.  the  radiant  transfer  of  heat. 

The  contribution  of  the  fourth  mechanism,  radiant  heat 
transfer,  was  assumed  negligible  in  all  the  papers  reviewed 
for  this  study.  In  addition,  the  variation  of  the  physical 
properties  with  temperature  and  pressure  was  neglected. 

The  Linear  Model 

Hadidi(6)  suggested  the  following  assumptions: 

1.  An  instantaneous  thermal  equilibrium  between  the  injected 
fluid  and  the  model  fluid  and  sand  is  always  attained, 
that  is,  the  convective  heat  transfer  coefficient  is 
infinite,  and 
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2a  the  heat  transferred  by  conduction  in  the  sand  and  in 
the  fluid  is  negligible. 

The  differential  equation  for  this  model  and  its  solution  are 


(PC  (1  -  6)  +  p  C 
r  pr  y  kw  pw 


9T  9T 

—  -  v  p  C  4>  — 

at  w  pw  ax 


(6) 


V, 


dx 

dt 


v  c  C  <t 
^w  pw 


nC  ( 1  d>  )  +  p  C  d> 
wr  pr  *  Mw  pw  y 


=  constant 


(7) 


XF  =  VFt 


(8) 


According  to  this  approach,  the  temperature  of  the  outlet  fluid 
would  vary  with  time  in  the  manner  illustrated  by  Curve  1  in 
Figure  3„  However  Preston  et  ai(ll)  found  that  the  correspond 
ing  outlet  fluid  temperature  profile  had  the  characteristic 
S-shape  shown  in  the  same  figure. 

Hausen (7)  and  Schumann (21)  independently  derived  equa¬ 
tions  describing  the  process  under  the  assumption  that  "mechan¬ 
ism  2",  heat  conduction  from  point  to  point  within  either  phase 
may  be  neglected.  The  corresponding  equations  are 


9  T  aT 

p  C  <p  — —  =  -  vp  C  d>  -  ha(T  -  T  ) 

Hw  pw  y  pw  V  w  r 

9  t 


(9) 


p  c  (l  - 

m  r  pr 


9T 


8t 


=  ha  (T  -  T  ) 
w  r 


10) 


The  solutions  to  Equations  9  and  10,  which  are  integrals  con¬ 
taining  Bessel  functions,  have  been  evaluated  by  numerical  or 
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graphical  methods.  Klinkenberg  (9)  found  that  the  exact  solu¬ 
tions  may  be  approximated,  within  fixed  limits,  by  probability 
integrals.  The  fit  was  made  by  choosing  the  appropriate  upper 
limits  of  integration.  Klinkenberg  proposed  the  approximations: 


T  -  T. 
w i 

Tn  -  T. 

1  l 


1  _  1  i 

1  +  erf  (/Z  +  -  -  /Y  +  - -)  !  (11) 


8/Z 


8/Y 


T  -  T, 
r_ l 

T,  -  T. 

1  l 


1  +  erf  ( /z 


__  -  /y  -  — -}  j  (12) 

8 /Z  8 /Y  ! 


where 


Y  = 


ha  x 


(13) 


C  p  v 
pw  w 


Z 


ha 


C  p  (1  -  <j>) 
pr  r 


(  t 


x 


v 


(14) 


=  water  injection  temperature,  CF 
Equation  12  may  be  used  for  any  value  of  Z,  Equation  11  for 
Z  >  1  provided  Y  >  2  in  both  cases.  The  maximum  error  under 

B 

these  conditions  is  less  than  -1%. 

Preston  et  al(19)  applied  the  Hausen  and  Schumann  model 
to  studies  of  packed  columns  to  determine  the  effect  of  the 
flow  rate,  the  column  length  and  the  sand  size  on  the  convective 
heat  transfer  coef ficient ,ha.  The  water  velocity,  based  on 
the  effective  flow  area,  varied  from  2  to  25  feet  per  hour,  the 
sand  size  from  20  to  170  U.S.  mesh  and  the  column  length  from 
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6.03  to  24  inches o  The  convective  heat  transfer  coefficient 
was  found  to  be  independent  of  the  sand  size  and  the  length 
of  the  column,  but  dependent  on  the  fluid  velocity  raised  to 
the  108  power,  according  to  the  equation 


ha  =  54.5  v1’82  (15) 

It  should  be  remarked  that  ha,  calculated  by  applying 
Equations  9  and  10  is  fictitious,  if  the  conduction  mechanism 
is  not  negligible. 

Jenkins  et  al(8)  essentially  made  the  opposite  assump¬ 
tion  concerning  the  mechanism  of  heat  transfer.  These  authors 
considered  conduction  to  be  an  important  mechanism  and  assumed 
negligible  heat  flow  by  convection  between  the  solid  and  the 
fluid.  Consequently,  the  system  can  be  described  by  one  dif¬ 
ferential  equation: 


3T 

(p  C  (1  -  <t>)  +  p  C  <j>)  — 
r  pr  Kw  pw 


3T  3  2T 

vp  C  4>  —  +  k  -  (16) 

KW  pWT  ,  .  o 

^  3x  3x^ 


where 


k  =  kw<j>  +  kr(l  -  <j>) 


The  solution  to  Equation  16  is 
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Equation  17  was  found  to  fit  the  experimental  results  of  Pres¬ 
ton  et  al  when  an  average  thermal  conductivity  of  1.0  BTU/hr. 
ft o °F  was  usedc  Any  contribution  of  the  convective  mechanism, 
which  may  be  present  in  this  value  for  the  thermal  conductivity, 
was  neglected. 

Kunii  et  al(10)  carried  out  experimental  work  on 
heat  transfer  in  linear  beds  of  glass  beads  and  sand  through 
which  fluids  were  flowing .  The  data  were  interpreted  in  terms 
of  an  apparent  or  effective  thermal  conductivity,  „  The  use 
of  an  effective  thermal  conductivity  implies  that  the  system 
possesses  a  single,  mean  temperature  at  every  pointo  The  mathe¬ 
matical  model  of  Kunii  et  al ,  however,  differs  from  that  of 
Jenkins  et  al  in  the  sense  that  kg  does  include  any  effect  of 
the  convective  heat  transfer  mechanism,  if  present.  Correla¬ 
tions  for  kg  vs  Reynolds  numbers  based  on  particle  diameter 
were  given  for  five  different  fluids:  helium,  nitrogen,  air, 
carbon  dioxide  and  liquid  water.  To  accurately  measure  values 
for  k  at  low  flow  rates  was  difficult.  This  is  specially 
true  for  the  sand-water  system,  for  which  only  one  correla¬ 
tion  for  modified  Reynolds  numbers  less  than  0.005  was  pre¬ 
sented.  The  values  of  kQ  were  found  to  increase  significantly 
with  the  mass  velocity  of  the  fluid. 

Green  et  al(5)  derived  a  general  set  of  differential 
equations  for  the  linear  model  with  all  three  heat  transfer 
mechanisms  included.  These  equations  are 
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For  the  water  phase 


3T 


P,  C  A 


w 


-  -  vp  C  <}> 


w  pw  '  w  pw 


3T  3ZT 

+  k  ♦  _ £ 
3x  W  3x2 


ha (T  -  T  )  (18 

w  r 


For  the  solid  phase 


3T 

prCpr(1  '  — 

r  pr  at 


3  ZT 


kr(l  -  <J>) 


3x 


— —  +  ha  (T  -  T  ) 
2  w  r 


(19) 


Equations  18  and  19  were  written  in  dimensionless  form  by 
choosing  the  parameters : 


a  = 


ha  ex 

( — ) 1/2  -s 


0) 


k  * 
w 


V 


T*  = 


T  -  T 


T,  -  T. 
1  r 


21) 


ha  , 

( - )  1/2  X 

V 


(22) 


^ia  1/2 

n  =  (— )1/z  vt 

k  (j> 
w 


(23: 


=  injection  temperature, 


and  solved  with  different  parametric  values  to  study  their 
effects.  The  systems  studied  were  water  flowing  through 
silica  beads  and  through  glass  beads.  At  A  values  of  0.114 
or  less,  the  numerical  solutions  were  found  to  approach  the 
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curves  assuming  =  kr  =  0,  implying  that  the  convective 
mechanism  essentially  controls  the  profile  shape.  When  X  was 
equal  to  0.342  or  higher,  the  numerical  solutions  more  closely 
approach  those  solutions  based  on  the  assumption  of  =  Tr. 

In  the  intermediate  X  range,  both  mechanisms  of  heat  transfer 
control  the  temperature  profile.  Green  et  al  have  explained 
how  the  above  numerical  calculations  may  be  used  to  select  a 
value  for  the  effective  thermal  conductivity,  k  ,  for  describ- 
ing  a  given  system  by  a  single  equation.  They  applied  the 
procedure  to  the  data  of  Preston  et  al(19)  and  found  that  kQ 
depended  on  the  fluid  velocity,  v.  A  relationship  between 
k  and  v  for  the  systems  silica-water,  silica-ethyl  alcohol 
and  glass-iso-octane  was  presented  in  graphical  form.  The 
velocity  range  studied  was  0  to  22  feet  per  hour.  The  great¬ 
est  deviation  between  calculated  and  predicted  k^ ' s  was  about 
20%. 

The  Cylindrical  Model 

The  multidimensional  physical  model  constructed  and 
operated  by  Malofeev (13 , 14)  had  a  sand  section  of  1  inch  high, 
28  inches  long  and  9  inches  wide.  The  surrounding  impermeable 
strata  were  simulated  by  clay  layers,  each  5  inches  thick. 

A  series  of  five  experiments  were  performed  at  different  hot 
water  injection  rates. 

Spillette (22)  calculated  the  thermal  behavior  of  Malo- 
feev's  physical  model  for  three  water  superficial  injection 
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velocities:  0.52,  2.07  and  4.14  feet  per  hour.  The  method 

of  the  characteristics,  as  applied  in  reference  17,  was  used 
to  solve  the  energy  balance  equation,  based  on  the  Jenkins 
et  al(8)  physical  model.  That  is,  the  convective  transport 
of  heat  was  accounted  for  by  the  assumption  of  instantaneous 
thermal  equilibrium  between  the  sand  and  the  fluid. 

This  literature  review  on  the  controlling  heat  trans¬ 
fer  mechanisms  in  porous  media  containing  flowing  liquids  has 
shown  that 

1.  there  is  a  difference  of  opinion  concerning  the  liquid 
velocity  at  which  the  convective  heat  transfer  mechanism 
essentially  controls  the  transport  of  heat,  and 

2.  the  concept  of  apparent  or  effective  thermal  conductivity 
provides  the  simplest  correlation  basis,  whether  the  con¬ 
vective  transport  of  heat  is  negligible  or  not. 
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THE  PHYSICAL  MODEL  FOR  BACK  PRODUCTION 


The  physical  model  was  represented  by  a  porous  cylinder 
of  radius  R  and  thickness  h  through  which  radial,  incompressible 
fluid  flow  occurred  in  a  horizontal  plane.  The  fluid  and  solid 
temperatures  were  assumed  equal  to  the  same  initial  constant 
value  throughout.  At  time  zero,  a  mixture  of  oil  and  water  at 
a  lower  temperature,  was  injected,  at  constant  pressure  and 
temperature,  into  every  point  of  the  external  cylindrical  sur¬ 
face.  The  following  assumptions  were  also  made: 

1.  The  impermeable  strata  had  zero  thermal  conductivity  in 
the  horizontal  direction  and  a  finite  thermal  conductivity 
in  the  vertical  direction. 

2.  The  permeable  reservoir  sand  had  a  finite  thermal  con¬ 
ductivity  in  the  horizontal  direction  and  an  infinite 
thermal  conductivity  in  the  vertical  direction. 

3.  The  solution  gas  in  the  reservoir  fluids  was  assumed  to 
be  negligible. 

4.  The  heat  transfer  process  within  the  permeable  sand  was 
described  by  using  an  average  effective  thermal  con¬ 
ductivity  selected  from  reference  5.  The  value  for  the 
effective  thermal  conductivity  was  held  constant. 

5.  The  fluids  entering  the  hot  cylinder  were  not  heated  by 
conductive  heat  transfer  through  the  sand  before  reach¬ 
ing  the  boundary  r  =  R. 

Values  of  the  physical  properties  were  estimated  at 
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average  conditions  of  temperature  and  pressure  and  were 
assumed  to  remain  constant.  The  oil  and  water  viscosities 
were  considered  to  be  known  functions  of  temperature. 

Other  assumptions  concerning  the  calculations  for  the 
heat  flux  through  the  flat  surfaces  of  the  cylinder  and  the 
fluid  distributions  are  given  in  the  following  pages. 

To  derive  the  energy  balance  equation,  consider  a 
cylindrical  element  of  cross-sectional  area  dA  and  height 
h,  which  are  depicted  schematically  in  Figures  4A  and  4C. 

The  total  energy  balance  may  be  written  as  follows : 


Heat  in  by 
conduction  + 
during  dt 


I-Ieat  in  by 
bulk  flow  * 
during  dt 


Heat  in  or 
out  by  conduction 
through  the  flat 
surfaces  during  dt 


Heat  out  by 
conduction  + 
during  dt 


Heat  out  by  Enthalpy  change  of 
bulk  flow  +  element  during  dt 
during  dt 


or  algebraically  in  the  form 


q  dt 


vpfcPf*2¥ 


(r  +  dr) 


hT  ,  ,  dt 
r-fdr 


+ 

qdt 


“  qr+drdt  +  vpfCpf  4>2  irrhT^dt  +  2-rrrhdr  d(pC^T) 
Rearranging 


=  2 wh  r  dr  d(pC  T) 

K  p 


±  qdt 


(24) 
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By  assuming  pC  =  constant  at  some  average  value  and  neglect- 

hr 

ing  the  term  containing  the  product,  dr  dt,  Equation  24  be¬ 
comes 

(qr  -  qr+dr>  +  vpfCpf*2„rh  (Tr+dr  -  Tr)  ±4 

3T 

=  27thpC  r  dr  —  (25) 

p  3t 

Expanding  the  expression  for  the  heat  flux: 

3T 

q  =  -  k  ( — )  2tt  rh  (26) 

r  e  3r 

in  a  Taylor  series  and  neglecting  all  but  the  first  two 
terms,  the  following  expression  for  qr+dr  is  obtained: 


dt  a  aT 

qr+dr  =  “  k  (— )  2Ttrh  +  —  [-  k  (-— )  27trh]  (27) 

r  e  3r  3r  e  ar 


Combining  equations  25,  26  and  27,  the  heat  balance  equation 
becomes : 


1  3  3T  3T  q  aT 

- [k  r(— )]  +  vpfC  —  +  -  =  pC  —  (28) 

r  )r  3r  p  3r  2whr  dr  p  3t 


By  choosing  an  average  constant  value  for  the  effective  ther¬ 
mal  conductivity,  k  ,  Equation  28  may  be  written  as 

a  2t  i  aT  •  aT 

a  -  +  (a  +M) - t  C  =  —  (29) 

ar2  r8r  at 
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C  = 

qT  - 


pC  h 
P 


q 

dA 


(31) 


(32) 


k 

—  (33) 


The 

following 

set 

of  boundary 

conditions 

were  chosen: 

1. 

T  = 

T. 

i 

=  constant 

when  r 

=  ^ 

2. 

3T 

0 

when  r 

=  0 

3. 

3r 

T  = 

T 

ST 

when  t 

=  0  for  0  .<  r  <  Rjj 

Calculation  of  the  C  Term 


To  calculate  the  heat  losses  or  gains  through  the  flat 
surfaces  of  the  cylindrical  model,  a  knowledge  of  the  tempera¬ 
ture  distribution  existing  in  the  overburden  and  in  the 
underburden  at  the  end  of  the  steam  injection  process  is 
necessary.  This  temperature  profile  is  a  function  of  R  and 
t,  implying  that  the  C  term  in  Equation  29  is  not  a  constant. 
The  impermeable  strata  was  divided  into  a  series  of  concentric 
cylinders  within  which  C  was  assumed  to  be  constant.  Equation 
29  was  solved  numerically  using  the  boundaries  of  these 
cylinders  as  grid  points.  This  approach  allowed  the  use  of 
different  C's  for  different  grids  and  thus,  took  into  account 
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the  dependence  of  C  on  R. 

The  computation  of  C  for  every  cylinder  was  done  by 
assuming  one-dimensional  heat  flow  from  the  sand  during  in¬ 
jection.  A  summary  of  the  relevant  equations  follows. 

Consider  a  cylindrical  element  with  an  inner  radius 
r j  and  an  outer  radius  rj+^  where,  j  and  j+1  are  subscripts 
indicating  the  j  and  the  (j+1)  grid  intersections  in  the  r- 
direction.  The  temperature  profile  in  the  cylinder  after 
injection  is  given  by  the  solution  of  the  equation 


3  2T  1  3T 

3Z2  a0B  3t 

with  the  boundary  conditions , 

1.  T  =  TgT  when  Z  =  0 

2 .  T  =  T .  when  Z  =  00 

1 

3.  T  =  T^  when  t  =  0  for  any  Z 

The  solution  is 


(34) 


T  -  T,  (T^m  -  T_.  )  erf c  ( - - 

1  ol  1  ~  / - r — 

2/aOBtj 


) 


(35) 


The  injection  time  for  each  cylinder,  t ^ ,  was  calculated  by 
the  method  of  Willman  et  al(23) ,  by  subtracting  the  time 
necessary  for  the  steam  to  reach  the  distance  r^  from  the 
combined  injection  and  soaking  time,  tj. 

Equation  34  was  then  solved  using  Equation  35  as  the 
initial  condition  to  determine  the  temperature  distribution 
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in  the  overburden  cylinder  during  the  back-flow  period. 
Hence,  the  new  boundary  conditions  for  Equation  34  are 


1. 

T  = 

S 

T 

S 

(t) 

when 

Z  =  0 

2. 

T  = 

T. 

l 

when 

Z  =  00 

3. 

T  = 

T. 

l 

+  (tst 

-  V 

erfc 

2/^OBEj 


)'  when  t  =  0  for  any  Z 


The  subscript  S  emphasizes  that  the  given  condition  is  valid 
at  Z  =  0  only. 

The  solution  is  given  by  Carslaw  and  Jaeger (3)  as 
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The  temperature  of  the  element  at  Z  =  0  (see  Figure  5) , 

0g (t)  =  Tg(t)  -  ,  must  be  expressed  as  a  function  of  time 

before  attempting  to  integrate  Equation  36.  0g (t)  was  approxi¬ 
mated  by  a  step  function  as  follows: 


es(t)  =  e(1) 

es(t)  =  e(2) 

es(t)  =  e(3) 


constant  for  0  <  t  <  x^ 
constant  for  <  t  < 

constant  for  T2  <  ^  <  T3 


es(t) 


=  0 


(N) 


=  constant  for 


XN-1< 


N 


Figure  6  shows  the  schematic  representation  of  the  step 
function  approximation.  The  second  integral  of  Equation  36 
becomes 
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The  first  integral  appearing  in  Equation  36  was  not  evaluated 
directly.  Instead,  its  derivative  was  calculated  using  inte¬ 
gration  by  parts.  The  algebraic  manipulation  of  the  equation 
led  to  the  following  expression  for  the  heat  flux,  q^,  (see 
development  in  Appendix  1) . 
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where  q^  was  taken  as  positive  if  the  heat  flux  was  in  the 
negative  (downward)  Z-direction. 

Assuming  that  the  underburden  behaves  in  the  same 

p- 

manner  as  the  overburden,  the  final  expression  for  qT  was 
found  to  be 
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where  t 


O  •  •  •  9 


total  back  production  time 

first  back  production  time  increment 

second  back  production  time  increment 

third  back  production  time  increment 

•  9 

th 

=  n  back  production  time  increment 


The  term  C  was  calculated  for  every  cylindrical  element  at 
every  time  step  using  Equations  31  and  39. 


Calculation  of  M 

The  term  M  appearing  in  Equation  29  was  defined  by 


Equation  30. 

In  field  units  it  becomes 
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Kr  pr  Y  Kf  pf 

(43) 

An  examination  of  Equations  40  to  43  showed  that 
expressions  for  the  rate  of  back  production  and  the  composi¬ 
tion  of  the  mixture  of  oil  and  water  flowing  had  to  be  deter¬ 
mined  before  solving  Equation  29  for  the  temperature  profile. 
Neither  of  these  two  variables  is  a  known  function  of  time 
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and  both  are  dependent  on  the  temperature  distribution  exist¬ 
ing  in  the  reservoir  at  any  time.  In  addition,  the  rate  of 
back  production  is  dependent  on  the  fluid  saturation  distribu¬ 
tion  prevailing  in  the  reservoir,  which  is  also  a  function  of 
time.  A  numerical  approach  was  chosen  for  solving  Equation 
29  which  facilitated  the  stepwise  change  of  values  of  these 
variables  with  time.  The  procedure  followed  will  be  explained 
in  a  step-by-step  fashion  for  the  sake  of  clarity, 
a)  Calculation  of  the  Initial  Rate  of  Back  Production 

A  schematic  representation  of  the  saturation  distribu¬ 
tion  during  steam  in jection ( 23)  is  shown  in  Figure  7.  The 
"actual"  saturation  distribution  may  be  approximated  by  taking 
average  values  in  each  zone  as  shown  in  the  same  Figure.  The 
initial  rate  of  back  production  was  calculated  from  the  equa¬ 
tion  : 
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This  equation  may  be  simplified  to 
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due  to  the  small  contribution  of  the  neglected  terms  for  small 
values  of  back  production 

b)  Composition  of  the  Fluids  Flowing  in  the  Hot  Cylindrical 

Reservoir 

Assuming  that  any  fluid  with  a  saturation  exceeding 
the  corresponding  irreducible  saturation  is  displaced  by  the 
fluids  flowing  immediately  behind,  the  composition  of  the 
fluids  flowing  at  every  zone  may  be  computed  from 
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The  fraction  of  water  or  oil  moving  in  the  hot  portion  of  the 
cylindrical  reservoir  was  obtained  by  averaging  the  corres¬ 
ponding  values  in  the  steam  and  the  hot  water  zone  "2"  (see 
Figure  7a)  on  a  volume  basis.  The  value  of  M  was  then  cal¬ 
culated  from  Equation  40  using  values  for  p  and  C  estimated 

r 

at  the  average  temperature  and  pressure, 
c)  Calculation  of  the  Temperature  Profile 

After  choosing  a  small  time  increment,  At,  for  which 
M  is  assumed  constant,  Equation  29  was  numerically  solved. 

The  temperature  distribution  was  calculated  in  the  following 
manner. 

The  hot  part  of  the  cylindrical  reservoir  was  divided 
in  m  elements  of  equal  (Ar) ,  as  shown  in  Figure  8b  (m+1  grid 
intersections)  and  expressions  for  an  equal  number  of  C's 
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were  established  according  to  Equations  31  and  39.  These 
expressions  are  of  the  form 

cp  =  a  -  b  T  'lj  (48) 

where 

a,  b  are  positive  constants 
(1) 

T  =  temperature  of  element  p  when  t  =  At 

P 

Equation  29  was  replaced  by  finite  differences  using  third 
order  correct  approximations  for  the  distance  derivatives 
and  the  Crank-Nicholson  approach  for  the  time  derivative. 

The  approximations  for  the  distance  derivatives  are 
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at  j  =  2  (See  Figure  8B) 
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and  Equation  29  becomes  at  r  =  0  (j  =  l) 

32T  3T 

( 2  a  +  M) -  =  —  (52) 

3r2  3 1 

Replacing  the  finite  differences  approximations  in  Equation 
29  and  applying  the  Crank-Nicholson  Method,  the  following 
expressions  were  obtained: 
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n  represents  the  time  level  at  which  the  temperatures 
and  coefficients  are  being  calculated.  For  the 
first  time  increment  n  =  0 
Applying  the  initial  condition. 
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By  writing  Equations  53,  56  and  57  (written  for  every  point 
from  j  =  3  to  j  =  m)  a  set  of  linear  algebraic  equations  was 
obtained.  In  matrix  notation,  the  set  of  equations  is  of  the 
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The  right  hand  side  of  this  system  is  known.  The 
elements  of  the  m  x  m  matrix  on  the  left  hand  side  may  be  cal¬ 
culated  using  the  appropriate  equations.  The  solution  for  the 
T- vector  yields  the  temperature  distribution, 
d)  New  Saturation  Distribution 

The  composition  of  the  fluids  flowing  in  the  reser¬ 
voir  is  essentially  controlled  by  the  saturation  distribution. 
In  addition,  since  the  cooling  rate  of  the  reservoir  sand  is 
dependent  on  the  composition  of  the  cooling  fluids ,  the  tem¬ 
perature  distribution  also  depends  on  the  saturation  profile. 
This  dependence  is  stressed  by  the  fact  that  the  heat  capacity 
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cf  oil  is  approximately  one  half  of  the  heat  capacity  of 
water. 

The  displacement  mechanism  for  mixtures  of  oil  and 
water  displacing  mixtures  of  oil  and  water  of  different  com¬ 
positions  at  different  temperatures  has  not  been  studied  yet. 
The  procedure  followed  in  the  calculation  of  the  saturation 
distribution  as  a  function  of  time  is  detailed  in  an  example 
problem  and  may  be  described  as  follows. 

The  hot  portion  of  the  reservoir  was  divided  in  m 
cylindrical  elements  of  equal  Ar  and  the  oil  (or  water) 
saturation  in  every  element  was  changed  in  a  stepwise  manner 
with  respect  to  time.  The  following  calculation  steps  were 
carried  out: 

1.  Values  for  the  irreducible  water  and  oil  saturations  were 
selected  from  known  relative  permeability  curves. 

2.  The  volume  of  the  liquids  with  saturations  in  excess  of 
the  pertinent  irreducible  saturation  was  calculated 
for  every  element  in  the  reservoir. 

3.  The  composition  of  the  fluids  being  displaced  from  every 
zone  was  calculated  using  Equation  46. 

4.  By  assuming  a  "piston-like  displacement"  mechanism,  the 
oil  saturation  in  the  last  element  of  the  steam  zone  was 
changed  to  a  higher  value 
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was  equal  to  the  total  displaceable  pore  volume  of  the  last 
element  in  the  steam  zone  and  that  of  the  hot  water  zone  at 
the  steam  temperature . 

Subsequent  stepwise  changes  in  the  oil  saturation  of 
the  steam  zone  were  made  in  a  similar  manner.  For  example, 
the  oil  saturation  of  the  second  element  was  set  equal  to  the 
saturation  given  by  Equation  61  when  the  total  volume  of  fluids 
produced  was  equal  to  the  total  displaceable  pore  volume  of 
the  first  and  second  elements  and  that  of  the  hot  water  zone 
at  the  steam  temperature.  A  piston-like  displacement  mechanism 
was  also  assumed  to  occur  in  the  cold  water  and  cold  oil  zones. 
The  oil  saturation  in  the  "f lushed-out"  steam  and  hot  water 
zones  was  subsequently  changed  to  that  of  the  cold  water  zone 
as  the  " f lood-f ront"  progressed  toward  the  wellbore.  Follow¬ 
ing  the  influx  of  the  cold  oil  into  the  coLd  water  zone, 
similar  changes  in  saturations  were  made  in  the  preceding  zones, 
e)  New  Production  Rate 

From  the  temperature  profile  obtained  at  the  previous 
time  step,  the  oil  and  water  viscosities  were  calculated  at 
the  average  temperature  of  every  element.  The  relative  per¬ 
meability  to  oil  and  the  relative  permeability  to  water  were 
calculated  for  every  element  using  the  saturation  distribution 
at  the  previous  time  step.  The  new  rate  was  obtained  from 
the  equation: 
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The  cold  water  zone  disappears  when  Rcq  =  .  Equation  62 

is  still  valid  provided  Rcq  is  set  equal  to  R^.  A  general 
flow  diagram  for  the  calculations  of  back  flow  is  given  as 
follows : 
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Discussion  of  the  Model 

M  3T 

The  term  ( - ) appearing  in  Equation  29  takes  into 

r  3r 

account  the  contribution  of  the  bulk-flow  heat  transfer 
mechanism  to  the  sand-cooling  process.  The  term  C  accounts 
for  the  non-adiabatic  behavior  of  the  flat  surfaces  of  the 
reservoir  model. 

Both  conductive  and  convective  heat  transfer  mechan¬ 
isms  are  incorporated in  the  analysis  by  using  the  effective 
thermal  conductivity  concept. 

However,  the  proposed  model  does  not  consider  multi¬ 
dimensional  heat  transfer  and  is  limited  to  constant  physical 
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properties  with  the  exception  of  fluid  viscosity  variation 
with  temperature. 

The  effects  of  neglecting  horizontal  heat  conduction 
in  the  impermeable  strata  during  the  injection  period  is  to 
increase  the  calculated  sand  temperatures  near  the  point  of 
injection  and  to  lower  the  calculated  radius  of  the  hot- 
reservoir  zone.  Avdonin (1)  and  Spillette ( 22)  have  compared 
the  effect  of  this  assumption  on  the  temperature  profile  and 
the  radius  of  advance  of  the  heat  front  for  hot  water  injec¬ 
tion  processes.  Their  mathematical  model , however ,  does  not 
consider  the  fluid  saturation  distribution  observed  by 
Willman  et  al(23)  in  laboratory  experiments  using  sand- 
packed  cores  and  linear  fluid  flow. 

The  assumption  of  zero  thermal  conductivity  within 
the  permeable  sand,  in  the  horizontal  direction,  during  in¬ 
jection,  will  have  effects  similar  to  those  just  mentioned. 

The  magnitude  of  these  effects  may  be  expected  to  be  less  pro¬ 
nounced. 

The  assumption  that  the  sand  and  the  fluid  may  be 
replaced  by  a  hypothetical  solid  of  thermal  conductivity,  k  , 
may  be  considered  as  an  accepted  procedure  to  describe  heat 
transfer  effects  in  packed  reactors.  The  calculation  of  this 
effective  thermal  conductivity  for  the  combination  of  fluids 
encountered  in  oil  reservoirs  offers,  however,  some  difficul¬ 
ties  due  to  the  lack  of  data  and  the  choice  of  a  proper  tech¬ 
nique  for  its  calculation. 
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Those  assumptions  concerned  with  the  saturation  dis¬ 
tribution  both  during  injection  and  back  production  may  affect 
the  results  seriously.  This  is  the  result  of  the  important 
role  which  the  composition  of  the  fluids  flowing  play  in  the 
back-production  calculations.  A  mathematical  formulation  to 
describe  the  saturation  distribution  and  hence  the  composition 
of  the  fluid  flowing  with  the  initial  temperature  and  satura¬ 
tion  distributions  of  the  present  model  would  be  complex.  It 
was  believed  that  this  formulation  was  not  worth  attempting 
in  the  present  case  unless  the  heat  transfer  simplifications 
made  in  establishing  the  physical  model  are  removed  or  proven 
to  be  valid. 
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EXAMPLE  CALCULATIONS 


To  illustrate  the  use  of  the  physical  model  for  the 
calculation  of  back  production,  two  example  problems  were 
run  on  a  7040  digital  computer.  The  computing  time  for  each 
example  was  approximately  15  minutes  for  a  total  production 
time  of  600  days  with  time  increments  of  2  days. 

The  effective  thermal  conductivity  for  the  oil-water 
sand  system  was  calculated  using  reference  5  as  follows: 
Average  readings  for  the  silica-water  and  silica-ethyl 
alcohol  were  obtained  from  Figure  9  for  v  *  2  ft. /hr.  and 

averaged  on  a  volume  basis  according  to  the  equation 


where 


X  k 


w  e 


w 


+ 


X  k 
o  e 


o 


effective  thermal  conductivity,  silica-water 


system 

k£  =  effective  thermal  conductivity,  silica-ethyl 
alcohol  system. 

Thus,  data  for  the  silica-ethyl  alcohol  system  were  assumed 

to  be  applicable  to  the  sand-oil  system  in  the  first  approxi 

mation  due  to  the  lack  of  the  proper  data.  The  variation 

of  k  with  the  fluid  composition,  X,  was  accounted  for  in 
^2 


the  calculations. 
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Example  1 

The  oil  viscosity  versus  temperature  curve  for  this 
example  has  been  plotted  in  Figure  10,  that  corresponding  to 
Example  2  appears  in  Figure  27.  The  viscosity  readings  at 
90°F  are  109  and  2300  cp  respectively. 

The  choice  of  viscosities  in  this  order  of  magnitude 
was  made  in  order  to  compare  the  effect  of  viscosity  on  the 
back  production  rate. 

The  initial  radii  and  oil  saturation  distribution  of 
the  back  production  period  are  shown  in  Figure  11,  permeabil¬ 
ity  data  are  plotted  in  Figure  9.  Additional  reservoir  and 
fluid  properties  are 
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(as  condensate) 

(p  )—  =  52.28  lbm/cu.ft. 
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=  58.45  x  0.882 

=  51.45  lbm/cu.ft. 
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(Cpw)jf 

1.03 

Btu/lbm°F 

^  Cpo  ^  T 

0.50 

Btu/lbm°F 

^Cpr^T 

0.23 

Btu/lbm°F 

1! 

CQ 

o 

1.09 

Btu/ft. °F  hr 

k  = 

eo 

1.20 

Btu/ft.°F  hr  (from  Ref.  5  for 

v  <  2  ft/hr. 

k  = 

ew 

1.70 

Btu/ft. °F  hr.  (from  Reference 

The  oil  viscosity  versus  temperature  data  were  fitted  by  the 
following  fifth  degree  polynomial 


V  =  10 

o 


where 


3  =  -0.65019612  +  0.44799295  x  10J(-) 

T 

3  1  7  1 

+  0.75749706  x  10  (-5-)  -  0.60803823  x  10  ( — =-) 

T  TJ 

9  ^  ii  1 

+  0.59698224  x  10  (~j)  -  0.18164653  x  10  (-=-) 

T  T5 

The  water  viscosity  versus  temperature  were  calculated  from 
the  equation ( 18) . 
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The  relative  permeability  data  presented  in  Figure  10  were 
fitted  by  sixth  degree  polynomials.  The  equations  are  as 


follows : 
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FIG. 9 

EXAMPLE  I. 

RELATIVE  PERMEABILITY  VS.  WATER  SATURATION 
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=  -0.06920184  +  0.12099786  x  101  x  S 
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Figures  12  to  21  show  the  step-wise  change  in  saturation  of 
the  reservoir  elements  as  a  function  of  the  total  fluid  pro¬ 
duction,  W  ,  for  subsequent  time  steps. 

P 

The  computer  program  is  given  in  Appendix  2. 


Results 

The  results  of  the  computer  programs  are  plotted  in 
Figures  22,  23  and  24.  Figure  22  shows  the  temperature  dis¬ 
tribution  as  a  function  of  the  time  of  back  production  in 
days,  the  center  line  at  the  left  corresponds  to  the  production 
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well  and  the  dashed  line  at  the  right  to  the  external  radius  of 
the  heated  reservoir  zone. 

The  results  for  the  rate  of  back  production  as  a 
function  of  time  are  given  in  Figure  23,  where  the  two  sets 
of  points  show  the  effect  of  changing  the  chosen  At  from  4 
days  to  2  days,  on  the  rate  of  back  production.  The  discussion 
of  the  above  results  appears  after  Example  2. 

The  results  of  Example  2  appear  in  the  same  order  as 
Example  1. 

Example  2 


The  initial  radii  and  saturation  distribution  are 
represented  schematically  in  Figure  25.  The  additional  data 
used  are 
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T.  =  74-F 

1 

-  P 

e  w 

1600  psi 

Total  injection  time  =  720 

<f>  = 

0.362 

(p  ) —  =  60.13  lbm/cu.ft. 

w  T 

Swc  ~ 

0.17 

( P  o>  T  =  50.32  lbm/cu.ft. 
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DISCUSSION  OF  RESULTS 


In  field  practice,  a  well  is  shut-in  to  allow  the 
steam  to  condense  and  the  pore  volume  of  the  steam  zone  to 
fill  with  liquid  before  the  back  production  cycle  is  initiated. 
Therefore,  the  calculation  of  the  external  radius  of  the  cold 
water  zone  was  based  on  the  assumption  that  the  content  of 
the  steam  zone  at  the  beginning  of  the  back  production  cycle 
was  condensate 0 

The  total  production  rate  before  steam  injection  for 
the  reservoir  conditions  of  Example  1  was  55  BPD .  The  initial 
rate  for  the  back-flow  period  was  found  to  be  384.3  BPD,  which 
represents  a  7-fold  increase  in  rate.  The  corresponding  in¬ 
crease  in  rate,  for  Example  2,  was  approximately  the  same. 

It  was  not  possible  to  compare  the  calculated  results  with 
any  field  results  since  most  of  the  reservoir  data  reported 
in  the  literature  were  incomplete.  Long (12)  reported  that 
the  lowest  ratio  of  improvement  of  the  initial  production 
rate  of  seven  wells  was  approximately  10,  whereas,  Boberg{2) 
reported  a  ratio  of  only  2.5. 

The  initial  rate  of  back-flow  essentially  depends  on 
the  radii  of  the  different  reservoir  zones  attained  during 
the  injection  period.  The  position  of  the  cold  water-cold  oil 
interface  also  plays  an  important  role  as  indicated  by 
Equation  45. 
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The  calculation  procedure  of  the  back  production  rate 
does  not  account  for  any  beneficial  effects  of  heating  on 
damaged  wellbores ( 2) .  In  addition,  the  assumptions  of  con¬ 
stant  sandface  pressure  and  radially-uniform  injection  of 
fluids  into  the  reservoir  were  made. 

The  scatter  of  points  in  Figure  23  is  the  result  of 
the  large  changes  in  the  oil  saturation  of  some  elements 
assumed  in  the  numerical  calculation.  In  this  example  the 
oil  mobility  is  more  dependent  on  the  relative  oil  permeabil¬ 
ity  than  in  Example  2.  To  illustrate  the  nature  of  this 
dependence,  consider  the  decrease  in  rate  when  t  -  24  days. 

The  change  in  oil  saturation  of  Element  7  from  0.49  to  0.69 
(see  Appendix  2)  corresponds  to  a  change  in  the  resistance  of 
this  element  to  fluid  flow  from  0.003  to  0.015.  This  sudden 
increase  in  resistance  causes  the  sharp  decrease  in  rate,  ob¬ 
served  in  Figure  23 0 

Two  distinct  regions  of  different  slope  may  be  clearly 
distinguished  in  Figure  23,  The  division  between  these  re¬ 
gions  occurs  at  approximately  t  =  150  days.  In  the  left-hand 
region,  the  back-flow  rate  decreases  drastically  with  time, 
while  the  decrease  in  the  right-hand  region  is  much  less  pro¬ 
nounced.  An  examination  of  Figure  24  reveals  that  t  =  150 
days  corresponds  to  the  time  when  the  production  of  water 
stops.  A  possible  explanation  of  this  discontinuity  in  the 
back  production  rate  versus  time  curve  is  the  change  in  the 
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cooling  rate  of  the  reservoir  due  to  the  change  in  saturation 
from  water  to  oil.  Since  specific  heat  of  oil  is  approximately 
one-half  of  that  of  water,  the  thermal  energy  transferred  to 
the  oil  is  reduced  to  one-half  of  the  previous  value. 

A  change  in  the  size  of  the  time  step  from  4  to  2  days 
produced  the  following  maximum  variation  of  results: 


a) 

Temperature 

3% 

b) 

Back-flow  rate 

3% 

c) 

Total  fluids  production 

2% 

where  the  results  for  At  =  2  days  were  used  as  the  basis  of 
comparison.  For  back  production  times  less  than  12  days,  how¬ 
ever,  a  temperature  variation  of  29%  was  encountered  at  grid 
points  near  to  the  external  boundary  of  the  model. 

It  was  also  observed  that  the  calculated  temperature 
of  the  first  grid  point  from  the  center,  independent  of  the 
size  of  the  time  step,  was  always  lower  than  the  correspond¬ 
ing  temperature  of  the  next  grid  point  during  the  first  three 
or  four  time  steps  (see  Appendix  2) . 

A  smoother  profile  would  be  obtained  by  incorporating 
the  term  C  into  equation  53. 

The  back  production  rate  curve  for  Example  2  did  not 
exhibit  the  scatter  of  points  encountered  in  the  same  curve 
for  Example  1.  The  reason  for  this  behavior  may  be  explained 
by  comparing  the  oil  mobilities  as  follows: 
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Example  1 
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Example  2 
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3000 
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2300  ( °F) 
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3000 
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2300 


u 
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Thus,  a  given  change  in  oil  saturation  changes  the  oil  mobility 
of  Example  1  approximately  4.3  times  the  corresponding  change 
in  Example  2.  This  means  that  the  error  introduced  in  Example 
1,  due  to  the  step  function  approximations  of  the  SQ  versus 
time  curve,  will  be  reflected  on  the  rate  versus  time  curve 
as  a  sharper  decrease  in  rate  than  the  corresponding  decrease 
in  rate  for  the  same  approximations  in  Example  2. 
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CONCLUSIONS 


The  examples  show  that  the  saturation  conditions  of 
the  reservoir  at  the  end  of  the  injection  period  essentially 
control  the  cooling  rate  of  the  sand  during  the  earlier  part 
of  the  back  flow.  During  this  part  of  the  cycle,  a  rapid 
decrease  of  the  back  production  rate  with  time  has  been  ob¬ 
served  in  field  operations ( 2)  as  well  as  in  this  study. 

Relatively  large  step-wise  approximations  for  the 
oil  saturation  in  the  reservoir  elements  will  produce  smooth 
rate-time  curves  for  those  cases  where  the  ratio  (k/yQ)  is 
close  to  unity.  The  size  of  these  approximations  becomes 
critical  for  very  favorable  mobilities. 
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RECOMMENDATIONS 

Work  in  the  following  areas  is  necessary  prior  to 
drawing  detailed  conclusions  on  the  proposed  model  for  evalu¬ 
ating  the  "Push-Pull"  method  of  steam  injection  as  well  as  on 
the  method  itself. 

1.  The  mathematical  treatment  of  the  injection  period  should 
be  revised  and  two-dimensional  heat  transfer  accounted 
for. 

2.  A  mathematical  approach  to  describe  the  displacement  of 
oil  from  the  sand  in  the  presence  of  thermal  effects 
should  be  developed. 

3.  A  method  for  predicting  the  effective  thermal  conductiv¬ 
ity  of  oil-water-sand  and  oil-water-natural  gas-sand 
systems  is  needed  to  describe  the  thermal  behavior  of 


the  reservoir. 
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NOMENCLATURE 


c,2 

area,  ft 

o 

differential  area,  ft 

combined  specific  heat  of  sand  and  fluids  at  con¬ 
stant  pressure,  Cp  =  Cpr(l  -  <j>)  +  C  f<|>,  Btu/lbm°F 

combined  specific  heat  of  oil  and  water  on 

volume-fraction  basis,  C  -  =  X  C  +  X  C  .  Btu/lbm°F 

pr  w  pw  o  po 

specific  heat  of  substance 
permeable  sand  thickness,  ft. 

convective  heat  transfer  coefficient,  based  on  a 
unit  volume  of  bulk  porous  media,  Btu/hr. ft^°F 

rate  of  heat  loss,  Btu/hr. 

rate  of  heat  injection  at  sandface,  Btu/hr. 
thermal  conductivity,  Btu/hr. ft. °F 
absolute  permeability,  md 

combined  heat  capacity  of  sand,  oil  and  water  as 
defined  by  Equation  4,  Btu/ft3°F 

number  of  time  increments  used  up  to  a  fixed 
production  time  t 

pressure,  psia 

pressure  difference  between  the  outer  and  inner 
radius  of  reservoir,  psi 

heat  flow  rate,  Btu/hr. 

heat  flow  rate  through  element  pf  area  dA  of  flat 
surfaces  of  cylindrical  model,  q  =  qTdA,  Btu/hr. 

total  heat  flux  through  flat  surfaces  of  cylindrical 
model,  Btu/hr. ft^ 

r-direction  coordinate,  ft. 

outer-well  radius,  ft. 
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‘w 


wellbore  radius,  ft. 


% 

S 

t 

t. 


t  . 
3 


T 

T 

Ts 

v 

VF 

W 

W 

P 

x 

X 

XF 

Y 

Z 

Z  ' 


•  *•  “ 

m 


"effective  radius"  as  defined  in  Reference.  1 

radii  of  elements  1,  2,  ...  m,  ft. 

fluid  saturation,  fraction 

back  production  time,  t  *  NAt,  hr. 

combined  injection  and  soaking  time,  hr. 

combined  injection  and  soaking  time  for  the  jth 
reservoir  element,  hr. 

temperature,  °F 

average  temperature,  °F 

surface  temperature  of  flat  surfaces  of  the  model,  °F 

fluid  velocity  based  on  "effective"  flow  area 

(A  =  cross-sectional  area/<J>)  ,  ft/hr. 
err 

velocity  of  "square"  heat  front,  ft/hr. 

back  production  rate  of  combined  oil  and  water, 
res.  bbls/day 

cumulative  fluids  production,  res.  bbls. 
x-direction  coordinate,  ft. 

volumetric  fraction  of  oil  or  water  flowing,  fraction 
x-coordinate  of  "square"  heat  front,  ft. 


dimensionless  distance,  Y  =  (ha/k  <f>) 

w 

Z-direction  coordinate,  ft. 
integration  variable 


1/2 


x 


Greek  Letters 


a 

A 


thermal  diffusivity, 
forward  difference 


=  k/pCp ,  ft  /hr. 


(ha/kw<|>)  1//2vt 


dimensionless  time. 
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X 

mobility  of  fluids,  * ^  =  (K  K/y  )  +  (K  K/y  ) 

md/c  1  co  ro  /Ho'co  v  rw  '  w'  co , 

y 

fluid  viscosity,  centipoise 

V  = 

dimensionless  group,  v  =  (2kQB/M' h^aQB) t^^2 

K 

integration  variable 

P 

combined  density  of  sand  and  fluids,  p=p  (1  -  <f> )  +  p -<j> , 
lbm/ft3  r  f 

Pf 

combined  density  of  fluids,  p  ~  =  X  p  =  X  p  , 
lbm/ft3  f  w  w  °  ° 

ii 

-H 

Q. 

density  of  i^  substance 

T  l'T2**  *  “ 

cumulative  back  production  time,  =  iAt,  hr. 

♦ 

effective  porosity,  fraction 

Subscripts 


C 

connate  saturation 

CO 

cold  oil  zone  of  reservoir 

CW 

cold  water  zone  of  reservoir 

e  = 

effective 

f 

fluid 

HWl 

hot  water  zone  assumed  at  initial  reservoir 
temperature 

HW2  =  hot  water  zone  assumed  at  steam  temperature 


• 

l  = 

initial 

t 

3 

jth  grid  intersection  in  redirection 

n  = 

time  increment  level 

o  =  oil 

oi  =  initial  conditions  of  oil  zone 

OB  =  overburden 

r  =  rock,  relative  or  residual 

ST  =  steam 

W  =  water 
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Overburden  and  Underburden  Temperatures 
Versus  Back  Production  Time 


Define 


e  = 


z=o 


P . D. E .  and  B . C. 


T 

T 


-  T. 

l 


ST 


-  Ti 


e  =  0  (t) 

S  S 


320 


3  Z- 


30 


K  3 1 
o 


(A) 


I.C.  t  =  0  0  =  0^  erfc  -—3 

2/otOBtj 


b.c.  t  >  o  eg  =  es  (t) 

e  =  o 


at  Z  =  0 
at  Z  =  °° 


Putting  0  =  U  +  V 

Equation  A  can  be  written  as  two  P.D.E. 


3  2  U 
3  Z  2 


3  U 


a0B  8t 
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°  2/^S 


B.C.  t  >  0 
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=  0 


at  Z  =  0 
at  Z  =  ~ 
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I.C.  t  =  0 


V  =  0 


B.C.  t  >  0 


v  =  eg(t) 


at  Z  =  0 


V  =  0 


at  Z  =  » 


The  solution  of  Equation  B  is 


(z  -  Z1)2 


D  = 


1  00  4a  t 

-  /  f(Z’)Me 

o 


-  e 


(Z  +  Z  1 ) 
4a0Bt 


}  dZ 


2/tt  a  __t 
OB 


(D) 


where 


f(Z')  =  8  erfc  - 
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The  solution  of  the  Equation  C  is  (Duhamel's  theorem) 


v  = 


2'/i'“ob  ° 


6S(X) 


(t  -  \) 


3/2 


dx 


Putting 


4a0B^t  ”  ^ 


,  the  above  equation  may  be 


written  as 
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The  above  expression  can  be  simplified  and  integrated  by 
parts  to  give 
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Consider  0  (t)  as  a  step  function; 
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Equation  E  can  be  integrated  now  to  give 


V(Z,t)  = 
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The  total  heat  flux  is  given  by 
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The  term  C  in  Equation  29  for  a  cylindrical  element  P  may  be 
written  as: 


C  =  C 
N  ^LHS 
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and  N  denotes  the  number  of  time  increments  chosen  for  a 
back  production  time,  t. 


Calculations  for  the  First  Time  Increment 


In  the  following  illustration  a  time  increment  of 
4  days  was  used.  No  appreciable  change  in  the  results  was 
found  when  this  time  increment  was  halved# 


a)  Initial  rate  of  back  production 
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(Figure  9) 
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cold 

water  zone 

S 

o 

O 

• 

cn 

VO 

k 

ro 

=  0.306 

(Figure 

Sw  = 

0.31  , 

k 

rw 

=  0.014 

T  = 

95°F 

y  =  0 . 

w 

7225 

y  =  100 

o 

(44) 
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cw 


0.306  x  612  0.014  x  612 

-  +  -  =  14.392 


100 


0.7225 


0.00708  x  25  x  400 


1  330  1  131.5 

-  In  -  +  -  in  - 

6.12  131.5  14.392  80.8 


*  384.3  BPD  (45) 

Irreducible  oil  and  water  saturations  (from  relative  permeabil 
ity  curves) : 


-  0.26 
or 


Swr  =  °*18 


(Figure  19) 


The  hot  cylinder  was  divided  in  ten  equal  Ar  elements.  The 
volume  of  oil  and  water  in  excess  fo  the  above  irreducible 
saturations  were  calculated  for  every  element,  these  volumes 
are  shown  in  Table  1.  The  composition  of  the  mixture  of  oil 
and  water  flowing  initially  in  every  zone,  calculated  accord¬ 
ing  to  Equation  46  are  also  shown  in  Table  1. 
b)  Composition  of  fluids  flowing  in  the  hot  portion  of  the 

reservoir. 

Using  Table  2  and  Figure  11, 
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Table  1 


Element  or  Zone 


Volume  of 
Displaceable 
Liquids 


Total  Volume 
of 

Displaceable 

Liquids 


Composition 

(X  ) 
w 


1  120 

2  357 

3  596 

4  837 

5  1072 

6  1320 

7  1545 

8  1776 

9  2015 

10  2040 

Hot  water  zone  1  2800 

Cold  water  22460 

Cold  oil  221750 


Table 

Injection  Time 


Element  (hr) 


1  720 

2  710 

3  690 

4  660 

5  620 

6  560 

7  485 

8  400 

9  285 

10  145 


120 

1.0 

477 

1.0 

1073 

1.0 

1910 

1.0 

2982 

1.0 

4302 

1.0 

5847 

1.0 

7623 

1.0 

9638 

1.0 

11678 

0.5 

14478 

0.5 

36938 

0.23 

258688 

0.0 

C  term 

assumed  negligible 

C0  =  0.336292  -  0 . 000904T9 ( 1 J 
2  z  (1) 
C3  =  0.334623  -  0.000904^  ' 

C4  =  0.331998  -  0.000904T4(1) 

Cc  =  0.328243  -  0 . 000904Tt- ( 1} 

3  D  ( l) 

Cc  =  0.321980  -  0 . 000904TC  K  1 

6  6  (1) 
Cy  =  0.312825  -  0 . 000904T7 

CQ  =  0.300033  -  0.000904TQ (1) 

8  °  (1) 
C9  =  0.276317  -  0.000904TgVJ-' 

C1Q  =  0.226922  -  0 . 000904T1Q C1) 
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c)  Calculation  of  the  temperature  profile 
1)  Computation  of  M 


pf  =  Xwpw  +  Vo  =  °*913  x  52.28  +  0.087  x  51.45 

=  52.208 


XwCpw  +  XoCpo  =  0,913  x  1-03  +  0-087  x  0.50 

0.9839 


prCpr(1  ’  +  PfCpf* 


140  x  0.23  (1  -  0.266) 
37.298 


+  52.208  x  0.9839  x  0.266 


M 


0.03720  x  384.3  x  52.208  x  0.9839 


25  x  37.298 


=  0.7875 


2)  Calculation  of  C 

The  hot  reservoir  and  the  adjacent  strata  were  divided 
in  ten  equal  Ar  increments.  The  injection  time  for 
every  element  was  calculated  a  explained  earlier 
using  Equation  1  of  Reference  23.  The  term  C  for 
the  first  time  increment,  obtained  by  combining 
Equations  31  and  39, is  given  by 


2k 


C  = 


OB 


hpC 


hpCp  /"aOB 


<T  -  T.  )  (—  -  — 

/At  /At  +  t 


T^1*  -  Ti  “I 


/It 


M  ?  o  noj  d&ijjqmoO 
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Table  2  shows  the  injection  time  for  every  element 
and  the  corresponding  equation  for  the  C  term. 
Equations  53  to  57  (the  latter  written  for  j  =  3  to 
j  =  10)  with  the  C  term  replaced  by  the  above  ex¬ 
pressions,  gave  a  system  of  simultaneous  algebraic 
equations.  The  system  of  equations  may  be  written 
in  the  compact  matrix  notation  as 

AT  =  r 

A  check  of  the  solution  was  made  by  multiplying  the 
found  T-vector  by  the  original  matrix  and  subtract¬ 
ing  the  new  independent  vector,  r' ,  from  the  original 
r.  The  maximum  deviation  was  found  to  be  less  than 
0.0001%. 

d)  New  saturation  distribution 

1)  Fluids  production 


W  =  W.  At  =  384.31  x  4  =  1537.24  bbls 

P  i 


2)  New  radius  of  the  cold  oil  zone 

-  5>61  x  1537.24 

R  =  VI 31,5^  -  _ 

co  it  x  0.266  x  25  (1  -  0.26  -  0.18) 

=  128.7  ft. 

3)  The  above  production  equals  the  volume  of  the  dis¬ 
placeable  liquids  from  the  section  of  element  9  lo¬ 
cated  inside  the  steam  zone.  The  saturation  of  this 
element  should  be  changed  to  (Equation  61) : 


0.20  +  (1  -  0.26)0.5 


0.47 
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The  saturation  in  the  other  part  of  the  hot  cylinder 
remains  unchanged  for  the  second  time  increment, 
e)  New  production  rate 

Following  the  method  outlined  in  the  preceding  pages,  the 
new  production  rate  to  use  in  the  second  time  increment 
for  the  calculation  of  the  temperature  profile  was  found 
to  be 

W  *  359.90  BPD 


Second  Time  Step 

A  new  value  for  M  is  computed,  now,  using  the  new 
rate  and  calculating  a  new  composition  for  the  fluids  flow¬ 
ing  from  the  new  saturation  distribution.  New  values  for  C 
should  be  calculated  according  to  the  equation 


2k 


C  = 


OB 


hpC 


h0Cp  /l,0OB 


(T  -  T.)  (— :  -  . ) 

-  bi  1  /2At  /  2  At  +  t.. 


ip(l)  _  ip  _  rp(2) 

—  ■  ■  ■■■—  +  .  . . — 

/2 At  /2At  -  At 

where  T^  is  the  unknown  temperature.  The  value  for  T^  is 
obtained  from  the  previous  temperature  profile.  Equation  29 
is  replaced  for  the  finite  differences  approximations  and  the 
new  system  of  algebraic  equations  solved  for  the  new  T-vector. 
The  rate  at  the  end  of  the  second  time  increment  is  found  using 
the  new  temperature  profile  and  the  new  saturation  distribution. 
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THE  COMPUTER  PROGRAM 
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Nomenclature  of  Computer  Input 


BCO 

=  A  __ 
CO 

mobility  of  fluid  in  cold  oil  zone,  md/Cp 

BCW 

=  X 

cw 

3= 

mobility  of  fluid  in  cold  water  zone,  md/Cp 

CONOB 

=  k  __ 
OB 

= 

overburden  thermal  conductivity,  Btu/ft.hr.°F 

CWO 

=  k 

ew 

=s 

effective  thermal  conductivity  sand-water 
system,  Btu/ft.hr.°F 

CPO 

=  Cpo 

ss 

specific  heat  of  oil  at  constant  pressure 
and  average  temperature,  Btu/lbm°F 

CP ROCK 

il 

o 

* 

ss 

specific  heat  of  sand,  Btu/lbm°F 

CPW 

ti 

o 

c 

= 

specific  heat  of  water  at  constant  pressure 
and  average  temperature,  Btu/lbm°F 

DSCO 

=  AS 

CO 

— 

1  -  S  -  S  =  total  saturation  change 

or  wr 

in  cold  oil  zone,  fraction 

DSHW 

ashw 

ss 

1  -  S  -  S  =  total  saturation  change 

or  wr  ^ 

in  cold  water  zone,  fraction 

DIFOB 

a  OB 

= 

2 

overburden  diffusivity,  ft  /hr. 

DPR 

=  AP 

ss 

pressure  drop,  psi 

H 

=  h 

s= 

permeable  sand  thickness,  ft. 

PERM 

«  K 

— 

absolute  permeability,  md 

POR 

=  <P 

ss 

porosity,  fraction 

RH 

ss 

"effective  radius" ,  ft. 

RHW 

ss 

radius  of  cold  water  zone,  ft. 

ROI 

•H 

O 

& 

II 

ss 

initial  radius  of  cold  oil  zone,  ft. 

ROO 

0 

Q. 

II 

— 

3 

oil  density  at  average  temperature,  lbm/ft 

ROROCK 

=  p 

r 

ss 

3 

sand  density,  lbm/ft 
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CO 
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W© 

owo 

r  P 

3* 

050 

*q 

at 

at 

wqo 

ss 

oo2* 

m 

se 

WH8  A 

WH8C3 

3K 

60  c 

acma 

«qa 

ss 

H 

X 

OS 

♦ 

SB 

SB 

HH 

WHH 

at 

SB 

XOOflOH 

II-3 


ROW 

— 

pw 

■ 

water  density  at  average  temperature,  lbm/ft 

RST 

ss 

rst 

ss 

steam  zone  radius,  ft. 

SO 

— 

s 

o 

=S 

oil  saturation,  fraction 

SOR 

ss 

Sor 

= 

residual  oil  saturation,  fraction 

SW 

ss 

Sw 

ss 

water  saturation,  fraction 

SWC 

ss 

Swc 

ss 

connate  water  saturation,  fraction 

TI 

S3 

Ti 

ss 

initial  reservoir  temperature,  °F 

TIMEIN 

ss 

ss 

combined  time  of  injection  and  soaking,  hr. 

TST 

ss 

T 

ST 

a 

steam  temperature,  °F 

XW 

a 

X 

w 

a 

volumetric  fraction  of  water  flowing 

XO 

ss 

X 

O 

= 

volumetric  fraction  of  oil  flowing 

Nomenclature  of  Computer  Output 

TEMPERATURE  PROFILE  =  temperature  at  each  grid  point  to 

count  from  the  wellbore,  °F 

SOLUTION  ERROR  =  relative  error  (percent)  made  in 

the  solution  of  the  matrix  for  grid 
point  in  question 

ELEMENT  BY  ELEMENT  RESISTANCE  =  value  of  the  term 

(1/A)  In  (Rm/Rm_1)  for 

every  element,  Cp/md 
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PLE  NO  1 

SOURCE  STATEMENT 


EUR f R AN  SOURCE  LIST 


SIBFTC  DECKA  NOL I  ST , NODECK 

ICO  FORMAT  (4F10.3/3F10.3/3F10.3/3F10.3/4F10.3/2FIO. 3/5F10.3) 

101  FORMAT  (1H  13F10.6) 

102  FORMAT  (  14X  *  6F1 0. 3/  1  AX  ,  6F  10 . 3/  14X  ,  6F  10  •  3/  14X 1 6F  1 0.  3/  14X  ,  3F  1 0.  3  ) 

103  FORMAT  (8F10.2) 

104  FORMAT  (  10F8.3) 

105  FORMAT  ( 1H  8F15.3) 

106  FORMAT  (2F10.3) 

107  FORMAT  ( 10F8.1) 

108  FORMAT! lHL13Xf29H INITIAL  BACK  PRODUCTION  RATE* * F 1 5 • 2  * 4H  BPD ) 

109  FORMAT! 14X * 32HI NCREMENTAL  AND  TOTAL  PRODUCT  ION, 2F 10* 1) 

110  FORMAT! 14X,23HSATURATI0N  01 STR IBUT I  ON/ 14X , 10F6. 3 ) 

1111  FORMAT ( 1HL13X*! 5HRE  SERVO I R  AF TER • F 1 0. I  * 5H  DAYS) 

602  FORMAT! 14X, 19H0RIGINAL  SATURAT ION/ 14X , 10F6. 3 ) 

603  FORMAT  ( 14X ,23H TOTAL  TIME  OF  INJECT  ION/ 14X*10F6.1) 

COMMON  M«  CLHS •  DT I  ME * R* PERM,RH* BCO * BCW *H , DPR  *  DEL TAR , T I , RR 1 , 

1  X(10)«Y(10)«B( 10)*S0(10)*T(15*150) »W(15Q)*NTIME*IXY*GIXY*HIXY 
DIMENSION  Cl(l0),C2(i0)*C12(l0),CCP(10*l50)*C(10*l50)*0C(10*l50)» 

1  PPP(10,l50),TIMEIN(10)*RTOTME(iO)»CP(10),CRHS(10),XCl(iO)*XC2(10) 
READ  ! 5 9  100 )  RS T , RO 1 * RHWt RH » SOR* DSHW* DSCO * ROW , ROO* ROROCK , 

1  CPW,CPO*CPROCK*DIFOB*CONOB*TST»TI *XW,XO*DPR,POR»H,BCO,BCW 
READ  (5*103)  PERMtSWC 
READ  (5*104)  (SO(I) *1=1*10) 

REAO  (5*106)  CWO*COO 

READ (5*107)  (TIMEIN( J)*J*1,10) 

READ  (5*106)  DT I ME 

WRITE (6*  102)  RS  T  *  RO I  * RHW* RH * SOR* DSHW* DSCO *ROW « ROO* ROROCK * CPW, CPO* 

1  CP ROCK  *  D I  FOB* CONOB  *TST*TI*XW«  XO*DPR*  POR  *H*  BCO*  BCW  *CWO*COO »  PERM 
WRITE  (6*602)  ( SO ( I ) *  I* 1  *  10 ) 

WRITE  (6*603)  (  TIME  IN  (  J  )  »  J=*l  *  10 ) 

V0LHZ*(3.1416*RST**2*P0R*H*( 1 . -SOR- SWC ) / 5 . 6 l ) 


CPF*XW»CPW+XO*CPO 

ROF*XW»ROW+XO*ROO 

R0CEPE*R0R0CK*CPR0CK*( l .-POR )+ROF*CPF*POR 
PR* (5.61*WP)/(3. 141 6*P0R*H*DSC0 ) 

R=SQRT (R0I**2— PR) 

IF  ( R« LE • RH )  RsRH 

WW*(0#00708*H*DPR)/( ( 1 . /BCO ) *ALOG ( 330. /R ) ♦ ( 1 . /BCW ) * ALOG ( R/RHW ) ) 
RR 1*0*  25 
T I ME*DT l ME 

DEN*SQRT(3.1416*DIF0B) 

4  RTIME*1./SQRT(TIME) 

CPT*2« • (CONOB/DEN )*RTIME 
DO  50  1*1*10 

RTOTME ( I )*l./SQRT(TIME*TIMEIN( I ) ) 

CP( I )  =  2. * (CONOB /DEN )*( (TST-TI ) • ( RT I ME-RTOTME ( I ) )*TI*RTIME) 

50  CRHS( I )*CP  (  I  )/( H#ROCEPE ) 

IF  (TIME.EQ.2.*DTIME)  GO  TO  5 
CLHS*CPT/(H«ROCEPE) 

WRITE  ( 6  *  108 ) WW 
W I  =  w  w 


WP*WI*DTIME/24. 

9  WF*(0.03720»WW)/(H»P0R) 
M=  1 C 


X  M  =  M 


12  1J  30HUO2  KAHIHOI 


I  OH 
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LE  NO  i  FURIRAN  SOURCE  LIST  DECKA 

SOURCE  STATEMENT 

CONG=XW*CWO+XO*COO 
DI AO=CUNO/ROCEPE 
XMM=*WF*ROF*CPF*POR/ROCEPE 
DEL  TAR=RH/XM 

Y  (  I  )=(2.*DI AO+XMM)/ (DEL  TAR* *2) 

X (  l  )  - Y ( 1 )+l./DT IME 
L  -M—  I 
XL  =  L 

DO  l  N3 1  ,  L 
XN  =  N 

RR-XN*DEL  TAR 
X(N+1)=DIAG/(DELTAR»*2) 

Y(N*1  )  =  (DIA0+XMM)/(6.*XN*0ELTAR**2) 

L  CONTINUE 

IF  (WW.NE.Wl)  GO  TO  10 
0WF*(0.03720»WI )/(H*POR) 

OXMM=OWF*ROF*CPF*POR/ROCEPE 

0C0N=C0NQ 

OD  I  A=sD  I  AO 

0XM=0DIA/(DELTAR**2 ) 

OYM=<ODIA+OXMM) / ( 6. * XL *DE L T AR* *2 ) 

B(1)=-TST*1./DTIME 
DO  2  J*2,L 

B( J)=-2.*TST/DTIME-CRHS< J) 

2  CONTINUE 

B ( M ) =  < 2.*0YM*0XM-(2./DTIME> ) *TST- ( OXM+2. *OYM ) * ( T ST  +  T I ) /2.- ( X <  M ) +2 . 
I* Y ( M ) ) *  T I-CRHSI M ) 

NTIME=1 
I  X  Y*  I 
G I  X Y  =  0 • 

H  I  XY  =  0  • 

RIXY=DTIME/24. 

WRITE  (6,1111)  RIXY 
CALL  TEMP 
CALL  RATE 

WP=WP+ ( W I +W ( 1 ) ) »DTIME/24./2* 

WW=W( IXY) 

11  ROF*XW*ROW+XO*ROO 
CPF=XW*CPW+XO*CPO 

ROCEPE=ROROCK»CPROCK*( 1 POR ) + ROF*C PF *POR 
PR=(5.61*WP)/<3.1416*P0R*H*DSC0) 

ENDCW=R0I**2~PR 
IF(ENDCW.LT.O.)GO  TO  143 
RESORT ( ENDCW ) 

IF ( R.GT.RHJGO  TO  44 
143  R=RH 

44  IF  (WW.LT.W(l))  GO  TO  12 
DO  3  J  =  2 ,  M 

XCK  J)=CP(  J)-CPT*T(  JfNTIME) 

C 1 C J)=XC1( J)/ROCEPE 

3  CONTINUE 
TIME=2.*DTIME 
GO  TO  4 

5  DO  6  J=2 , M 

XC2 ( J  )  =CP ( J  )  f CP  T * ( SORT (2. )- 1 . ) *T ( J , NT  I  ME ) /SORT  (  2  .  ) 
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le  NO  1 

SOURCE  STATEMENT 


T  OK  1  RAN  SOURCE  LIST  DECK  A 


6  C2 ( J  )  =XC2 ( J )  /ROCEPE 
DO  7  J  =  2  *  M 

C 1 2  C J)=Cl( J)*C2( J) 

7  CONTINUE 
CLHS^CPT/RUCEPE 

B(  1 )=-Y( l)»T(2«NTIME)+(Y( I ) -1 . /DT I  ME ) *T ( I, NT  I  ME) 
B(2)=-(X(2)-2.*Y(2))*T(19NTIME ) - ( 2 . /DT I ME-2 . *  X ( 2 ) -3 . * Y ( 2) )* 

1  T (2*NTIME)~(X(2)+6.*Y(2)  )*T(3»NTIME)+Y(2)*T(4tNTIME)-C12(2) 

DO  8  J*3,M 

B(J)*-Y(J)*T(J-2fNTIME)-(X(J)-6.*Y(J))*T(J-ltNTIME)-(2./DTIM£- 
l  2.*X( J)+3.*Y(J) )*T( J,NTIME)-(X( J)+2.*Y( J) )*T(J+1»NTIME)-C12( J) 

8  CONTINUE 
GO  TO  9 


10  0(M)*B(M)— (X(M)+2.*Y(M)  ) *  T ( M+  l  *  NT  I  ME ) 

NT  I ME  =  NT I  ME  +  1 
I  XY=*  I X  Y+  l 
F I X Y- 1 XY 

RIXY=DTIME*FIXY/24. 

IF  (HIXY.EQ.GIXY)  WRITE  (6»liil)RIXY 
IF  (IXY.EQ.150)  GO  TO  32 
CALL  TEMP 

IF  IT! 7,NTIME).LE.TI )  GO  TO  32 
IF  (ABS(T(8,NTIME)-TI).LE.2.)  GO  TO  39 
IF  (ABS(T(9,NTIME)-TI ).LE.2.)  GO  TO  40 
IF  (ABS(T( IO,NTIME)-TI I.LE.2. )  GO  TO  41 
43  CALL  RATE 

DWP=( WINTIME- 1) +W(NT I  ME) \ *DT I  ME/ ( 2. *24. ) 
WP-WP+DWP 

IF  (HIXY.EQ.GIXY)  WRITE  ( 6 # 109 ) DWP , WP 
IF  ((WP. GE. 1535.). AND. (WP.LT. 3311.))  GO  TO  77 
IF  ( (WP.GE. 3311 .). AND. (WP.LT. 4856. ) )  GO  TO  33 
IF  ( (WP.GE. 4856. ). AND. (WP.LT. 6176. ) )  GO  TO  24 
IF  ( (WP.GE. 6176. ) .AND. ( WP.LT. 8085. ) )  GO  TO  28 
IF  ( (WP.GE. 8085. ) .AND. ( WP.LT. 8562. ) )  GO  TO  29 
IF  ( (WP.GE. 8562. ) .AND. (WP.LT. 9158.) )S0(3)=.49 
IF  ( (WP.GE. 9158. ). AND. (WP.LT. 9515. ) )  S0(2)=.49 
IF  ( (WP.GE. 9515. ) .AND. (WP.LT. 10177. ) )  GO  TO  30 
IF  ((WP.GE. 10177.). AND. (WP.LT. 11497.))  GO  TO  34 
IF  ( (WP.GE. 11497. ). AND. (WP.LT. 12569. ) )  GO  TO  35 
IF  ((WP. GE. 12569.). AND. (WP.LT. 14002. ) )  GO  TO  36 
IF  ( (WP.GE. 14002.  )  .AND. (WP.LT.  14479. ) )  GO  TO  205 
IF  ( (WP.GE. 14479. ). AND. (WP.LT. 14836. ) )  S0(2)=.69 
IF  ( (WP.GE.  14836. ). AND. (WP.LT.  16337. ) )  GO  TO  206 
IF  ( (WP.GE.  16337. ). AND. (WP.LT. 27316. )  )  GO  TO  207 
IF  ( (WP.GE. 27316. ). AND. (WP.LT. 31 107. ) )  GO  TO  114 
IF  ( (WP.GE. 31107. ). AND. (WP.LT. 33972. ) )  GO  TO  300 
IF  ( (WP.GE. 33972. ). AND. (WP.LT. 35881. ) )  GO  TO  302 
IF  ( (WP.GE. 3588 l. ). AND. (WP.LT. 36884. ) )  GO  10  304 
IF  ( WP.GE. 36834. )  GO  TO  307 
27  WW  =  W (  I  X Y ) 


X0=l.-XW 

IF  (HIXY.EQ.GIXY)  WRITE  (6»110)  ( SO l I ) » I = 1 , M ) 
GO  TO  ll 

12  CLHS=CPT/ROCEPE 


LL=NTIME*1 
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Lfc  NO  1  f  OKI  KAN  SOURCE  LI  SI  DECKA 

SOURCE  STATEMENT 

00  22  JJ=NTIME,LL 
X J J= J J 
DO  21  J  =  2  »  M 

CCP  (  J,  JJ  )=CLHS*  I  (TST-T I ) *  I  1  ./SCR  T  (  X  JJ  J-SQRI  (  DT  I  ME  )  /  SUR  I  (  X  J  J  *0  l  I  ME  ♦ 
i  TIMEINT J) ) )*Tl /SORT (XJJ) ) 

CC  =  0. 

DO  18  N  =  2 , NT  I  ME 

K=LL-N 

XNN  =  N 

CC  =  CC  +  CLElS*(SQRTlXNN>-SQRTt  XNN-1. ) )  * T I J , K  ) /SORT  I XNN* I XNN- 1 . ) ) 

10  CONTINUE 

21  C< J,JJ)=CCPI J,JJ)+CC 

22  CONTINUE 

DO  23  1=2, M 

PPP( I,NTIME)  =  C(  I, NT  I  ME  )-CLHS*T 1 1 , NT  I  ME ) 

2  3  DC  I  I , NT  I  ME ) =PPP ( I , NT  I ME ) +C (I *L L ) 

2  5  B< 1  )=-Y(l)*T(2tNTIME)  +  (Y(l)-l./DTIME)*T( 1, NT  I  ME) 

B(2  )=-IX( 2 )-2.*YC2)l*T( 1, NT  I ME >- 1 2. /D T I ME-2. *X I  2 > -3 . *Y I  2 ) )* 

1  T (2,NTIME)— (X( 2  )  +6  •  *  Y  I  2  )  )  *T ( 3 , NT IME ) ♦ Y ( 2 ) * T I 4 , NT  I  ME ) -DC I 2 , NT  I  ME ) 
DO  26  J=3 , M 

Bt J)=-YI J)*TI J-2,NTIME)-IXI J)-6.*Yt J)  ) *T ( J- 1 , NT  I  ME ) - ( 2 . /DT I  ME- 

1  2.*XIJ)+3.*YIJ) )«T(J,NTIME)~(X(J)+2.*Y(J) )*TIJ*1,NTIME)- 

2  DC ( J , NT  I ME ) 


26 

CONTINUE 

GO 

TO  9 

77 

SO ( 9 ) = .49 

X  w  = 

.841 

GO 

TO  27 

33 

SO  I  8 )  =  • 49 

SO  I  9 )  =  .49 

X  W  = 

.78 

GO 

TO  27 

24 

DO 

79  1=7, 

9 

79 

SOI  I  )  =  • 49 

SOI  10)  =  . 69 

x  w= 

.671 

GO 

TO  27 

28 

DO 

80  1=6, 

8 

80 

SOI  I )  =  . 49 

SO  I  9 )  =  • 69 

X  w  = 

.627 

GO 

TO  27 

29 

DO 

81  1=4, 

7 

81 

SOI  I  )  =  • 49 

DO 

82  1=8, 

10 

82 

SOI  I )  =  • 69 

X  w  = 

.44 

GO 

TO  27 

30 

DO 

83  1=1, 

7 

83 

SOI  I  )  = . 49 

DO 

84  1=8, 

10 

84 

SOI  I  )  =  .69 

XW  = 

.404 

GO 

TO  27 

34 

DO 

112  1=1 
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112 

113 


35 

115 

116 


36 

117 

118 


205 

120 


206 

122 


207 


114 

119 


300 

301 


302 

303 


304 

305 


30  7 
308 


39 

201 


40 

202 


l  I  (JlU  KAN  SOURCE  LIS 

SOURCE  STATEMENT 


SU(  I  )  =  •  4  9 
00  113  1=7,10 
SOI  I  )  =  .69 
Xto=  •  3  7  1 
GO  TO  27 
00  115  1=1,5 
SOI  I  )  =  .49 
00  116  1=6,10 
SOI  I )  =  . 69 
X W= • 326 
GO  TO  27 
DO  117  1=1,4 
SOI  I  )  =  • 49 
DO  118  1=5,10 
SOI  I )  =  . 69 
X W= • 29  3 
GO  TO  27 
00  120  1=3,10 
SO  I  I )  =  . 69 
X W=  •  247 
GO  TO  27 
DO  122  1=1,10 
SOI  I  )  =  . 69 
X W= • 2  1 9 
GO  TO  27 
SOI  10 1  =  • 82 
XW= • 209 
GO  TO  27 
DO  119  1=8,10 
SOI  I  )  =  .82 
X W= •  1 1 6 
GO  TO  27 
DO  301  1=6,10 
SOI  I  )  =  .82 
XW= • 065 
GO  TO  27 
DO  303  1=4,10 
SOI  I  )  =  .82 
XW= . 02 1 
GO  TO  27 
DO  305  1=2,10 
SOI  1 )  =  • 82 
XW=0.0 
GO  TO  27 
00  308  1=1,10 
SOI  I  )  =  .82 
XW=0.0 
GO  TO  27 
00  201  1=8,10 
T I  I , N  T I  ME )  =  T I 
M  =  7 

GO  TO  43 
DO  202  1=9,10 
T I  I ,NT I  ME )  =  T  I 
M  =  8 
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WX 
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MX 

.. 

ss.  i  me  ioe 
edo.=wx 
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VS  OT  00 
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LE  NO  l  HJKJRAN  SOURCE  LIST  DECK  A 

SOURCE  STATEMENT 

GO  TO  43 

41  T(lO»NTIME)=TI 
M  =  9 

GO  TO  43 
32  CALL  EXIT 
END 


lyHMSIATe  3J>IU02 


1  OM  3J 


!  i  *  l  >  i  •  1  Jr 

l*  OT  03 
11X3  JJA3  S£ 


' 


FORT KAN  SOURCE  LIST 


LE  NO  1 

SOURCE  STATEMENT 

SIBFTC  SOLUT  NUL I  ST , NODECK 
SUBROUTINE  TEMP 

199  FORMAT  ( IH  13F10.6) 

200  FORMAT  ( 14X, 19H TEMPERATURE  PROF ILE/ 14X , 1 1F6. 1 ) 

201  FORMAT  (  14X  »  1 4H SOLUTION  E RROR/ 14X , 1  OF  6 . 3 ) 

COMMON  M,CLHS,  DT  I  ME  ,  R  ,  PERM ,  RH  ,  GCO  ,  BCW  ,  FI ,  DPR  ,  DEL  T  AR  ,  T  I  ,  RR  1 

1  X(  10) *Y( 10) ,B(  10) ,  SO  (  10) *T ( IS, 150)  » W  (  150) *NTIMEf IXY 
DIMENSION  ERROR ( 10 ) 

DOUBLE  PRECISION  A ( 10, 10 ) , AA t 10, 10 ) ,BB ( 10 ) , XQT ( 10 ) , TE ( 1 1 , l 50 ) , 

1  D I F ( 10 ) , PCENT ( 10) ,AL( 10, 10) , AM( 10, 10) ,AU( 10, 10) ,AV( 10, 10) , XIS 

2  VALUE, 2X,2XX,ZV,ZW,ZQ,ZR, BOB 
DO  10  J=l,10 

DO  10  1=1,10 
10  A  I  1 , J)=0. 

A ( 1  ,  1 ) =-X ( 1 ) 

A(2,1)=X(2)-2.*Y(2) 

A ( 3 , 1 ) =Y ( 3 ) 

A ( 1 , 2 ) =Y (  1) 

A(2,2)=-2./DTIME~2.*X(2)-3.*Y(2)-CLHS 

A(3,2)=X(3)-6.*Y(3) 

A ( 4 , 2 ) =Y ( 4 ) 

A(2,3)=X(2)+6.*Y(2) 

A  I  3,3)=-2./DTIME-2,*X< 3)+3.*Y( 3)-CLHS 
A ( 4 , 3 )  =  X I  4 )— 6«*  Y ( 4 ) 

A ( 5 , 3 ) =Y ( 5 ) 

A ( 2 , 4 ) =—Y ( 2 ) 

A( 3,4)=XC  3)+2.*Y(3) 

A(4,4)=-2./DTIME-2.*X<4 ) +3. »Y < 4 )-CLHS 
A(5,4)=X(5)-6.*Y(5) 

A ( 6 , 4 ) =Y ( 6 ) 

LL=M-3 

DO  5  1=4, LL 

A (  I  ,  1  +  1  )=X< l)+2.*Y(  I ) 

A ( I  +  i, I+l)=-2./DTIME-2.*X( I  ♦ 1 ) +3. *Y U ♦ 1 ) -CLHS 
A( 142, I+l )=X( I+2)-6.«YC 1+2) 

AII+3,  1  +  1  )=Y<  1  +  3) 

5  CONTINUE 

A { M-2 , M- 1 ) =X ( M- 2 ) +2 . *Y ( M-2 ) 

A ( M— 1 , M— 1 )  =— 2 • / DT I ME-2 • *X ( M- 1 ) +3 • * Y ( M- 1 ) — CLHS 
A(M,M— 1)  =  X(M )— 6#  *Y( M ) 

A  (  M—  l  ,  M  )  =  X  (  M—  1 )  ♦  2  •  •  Y  I  M-  1 ) 

A(M,M)=-2./DTIME-2.*X<M)+3.*Y(M)-CLHS 
DO  7  J= 1 , M 
DO  7  1=1, M 

7  A A ( I , J )  =  A ( I,J) 

DO  8  1=1 ,M 

8  BBII )=B( I ) 

00  13  I = 1 , M 

13  XOT (  I  )  =  B (  I  ) 

C  CALCULATION  OF  L  AND  U 

DO  105  K  =  2  ,  M 
KK=K- 1 

DO  104  J  = 1 , KK 
AL ( J,K )=0.0 
AM(J,K)=0.0 


•v.;.  * 
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>LE  NO  l  FORTRAN  SOURCE  LIST  SOLUT 

SOURCE  STATEMENT 

AUIK, J )=C.O 

104  AV(K,J)=0.0 

105  CONTINUE 

IF ( A( 1 « 1 ) .EQ.O. )  GO  TO  150 

106  X  I $*A8S ( A (  1,1)) 

AL(  i,i)=SQRT(XIS) 

AU(  1 , l )  =  A ( 1,1)/ AL( 1, l) 

DO  107  J=2 , M 

AU ( l , J ) = A ( l,Ji/AL(l,l) 

107  AL(J9l)-A(J,l)/AU(i,l) 

DO  115  K=2 , M 

KK=K- 1 
VALUE=0. 0 
00  108  J-1,KK 

108  VALU£=VALUE+AL(K,J)*AU( J,K) 

ZX=A(K,K)— VALUE 

I  F ( ZX . EQ*  0 .  )  GO  TO  150 

109  ZXX=ABS ( ZX ) 

AL(K,K)=SQRT(ZXX) 

AU(K,K)=ZX/AL(K,K) 

KP-K-f  1 

KK=  K- 1 

IFIKP.GT.M  )  GO  TO  115 

DO  112  I -KP , M 

ZV=0.0 

ZWx0.0 

DO  110  LP= 1 , KK 
ZV=ZV+AL(K,LP)*AU(LP,I ) 

110  ZWXZW+AL(  I  ,LP)*AU(LP,K) 

AU ( K, I )=( AIK, I)-ZV)/AL(K,K) 

112  AL< l,K)xI A( I,K)-ZW)/AU(K,K) 

115  CONTINUE 

C  L  AND  U  ARE  CALCULATED 

C  PROCEEDING  TO  CALCULATE  L-INVERSE  AND  U-INVERSE 

DO  119  K=  l ,  M 
AMIK,K)-1*0/AL( K»K) 

119  AV(K,K)*1.0/AU( K,K) 

DO  125  K-2 , M 
KK=K-1 

DO  122  J  =  1 ,  KK 
ZC=0*0 

DO  120  La J  ,  KK 

120  ZQXZQ+AL(K,L)*AM(L,J ) 

122  AMI K,J )=-ZQ/AL( K,K) 

125  CONTINUE 
I JN-M  + l 

DO  135  KLX l , I JN 
K  -  I  JN-KL 

IFIK.LE.l)  GO  TO  136 

126  DO  130  JK= 1 , K 
J=K— JK 

IF(J«LT«1)  GO  TO  135 

127  ZR-0 • 0 
JP=J*i 

DO  128  L= JP ,K 
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>LE  NO  1  FOR THAN  SOURCE  LIST  SOLUT 

SOURCE  STATEMENT 

128  ZR=ZK*AU( J»L)*AV(L,K ) 

130  AVI J,K )=-ZR/AU( J,  J) 

135  CONTINUE 

C  PROCEEDING  TO  CALCULATE  G-  I  NVERSE-U- I NVERSE*L- INVERSE 

136  DO  140  K=  l » M 
DO  140  J=1,M 
A( J,K)=0.0 
DO  137  L=  l » M 

137  A(J,K)=A(J,K)+AV(J,L)*AM(L,K) 

140  CONTINUE 

GO  TO  152 

150  WR I TE ( 6# 151 ) 

151  FORMAT ( 1H  ,29H  THE  METHOD  IS  NOT  APPLICABLE) 

152  CONTINUE 

DO  153  1*1, M 
TE ( I • NT IME )=0« 

DO  153  J-  1 » M 

15  3  TE (  I,NTIME)=TE( I, NT  I  ME ) ♦A ( I ,J)*XQT( J) 

DO  14  1=1, I JN 

14  T ( I ,NT IME ) =SNGL (TE(I,NTIME) ) 

T ( M*l , NT IME )  =  Tl 

WRITE  (6,200)  ( T( I , NT  I  ME ) , I  =  1 , 1 JN ) 

DO  9  J= 1 , M 
B0B=0« 

DO  12  1=1, M 

12  BOB=BOB+AA( J,I)*TE( I,NTIME) 

DIF(J)=BB( JJ-BOB 

PCENT( J)=DIF( J)#100./BB( J) 

9  CONTINUE 
DO  15  J= 1 , M 

15  ERROR ( J ) =SNGL ( PCENT ( J ) ) 

WRITE  (6,201)  ( ERROR ( J  )  , J= 1 , M ) 

202  RETURN 
END 
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LE  NO  l 

SOURCE  STATEMENT 


FORTRAN  SOURCE  LIST 


$  I  BE  TC  NINAL 
SUBROUT 
200  FORMAT 
5000  FORMAT 
500  L  FORMAT 

5002  FORMAT 

5003  FORMAT 

5004  FORMAT ( 

5005  FORMAT ( 
COMMON 

1  X(10), 
01  MENS  I 
PWMAX= ( 


NOLIST 
INE  RATE 
(1H  L3FI 
(  IX,  13F6 
( 2X#  18F4 
(  l X , 7F  10 
(  IX, IF5. 
14X, 4HRA 
14X,29HE 
M,CLHS, 
Y(  10)  , B( 
ON  T  A  I  20 
1 . / BCO ) * 


, NODECK 

0.6) 

.1/5F6.  1  ) 

.2) 

.3) 

3, 12F6.2/5F6.2) 

TE ,F 10. I  ) 

LEMENT  BY  ELEMENT  RE S I S T ANCE/ 14X , 10F6 . 3 ) 

DTIME,R,PERM,RH, BCO,BCW,H, DPR, DEL  TAR, TI ,RR1 
10),S0( 10) ,T( 15,150) ,W( 150), NT  I  ME, I X Y ,G I XY , H I  X Y 
)  ,SW(20) tVISCOl (20) ,VISW(20) ,BCHW(20) , BBC (20) 
ALOG( 330./R)+( 1 . /BCW  )  «  ALOG  (  R/RFJ ) 


t 


XM=M 


T ( M  + 1 , NT  I  ME )  =  T I 
DO  1  J= 1 , M 
I  J  =  M  +  l- J 

SW(  I  J  )  =  l .-SO ( I J ) 

XKR0=-0. 69201 83 5E-0 1+0. 12099786E+0i*S0 ( I J )-0. 8639229 IE +  01* SO ( I J )*« 

1  2+0.28682803E+02*S0( I J )  * *3-0. 43337840E  +  02*S0 ( I J ) **4+0 . 30396 180E+ 

2  02*S0 ( IJ)**5-0.64854211E+0i*S0(IJ)**6 

XKRW=-0. 4883480 3E-0 1+0. 78 1 1 338 5E+00*S W ( I J )-0. 46846743E+0 1 *SW ( I J )** 

1  2+0.1 276 6994E+02*SW( IJ)*»3-0. 1289301 7E+02* SW ( I J ) **4+0 . 3556431 0E+ 

2  Ol*SW( IJ)**5+0. 150 0614 0E+0i*SW(IJ)**6 
TA(  I  J  )  =  (  T  ( IJ  +  1,NTIME)  +  T(IJ,NTIME) )/2. 

P0LYV0=-0.65019612E+00+0.44799295E+03*l./TA( I J ) +0. 75749706E+03* 

1  l./TA( I J ) **2-0 .608 0382 3E +07* 1 ./TA ( I J ) **3+0. 59698224E+09* 1 . /T A ( I J ] 

2  **4-0. 18 164653 E+li* l./TA( I J)**5 
V  I  SCO 1 ( I J ) -EXP(  2.30258509«P0LYV0) 

503  XM0B0=XKR0*PERM/VISC01 ( IJ ) 

VISWA=.55556*TA( I J) -26. 2 128 
VISWB-VISWA**2 
VISWC=8078.4+VI SWB 
VISWD*VISWA+SQRT( VISWC ) 

DENVIS=2.1482*VISWD-120. 

VISW( IJ)^100./DENVIS 
XMOBW-XKRW*PERM/VISW( I J) 

BCHW( I J )  =XMOBO+  XMOBW 
BBC (  I  J  ) s 1 • /BCHW ( I J) 

IF  (BBC( IJ).LE. 0.001 )G0  TO  504 
l  CONTINUE 


504 

1  x  = 

I  J-  l 

IF 

(  IX. EQ 

.0 

)  GO 

TO  506 

DO 

77  K  =  1 

.IX 

77 

BBC 

( K ) =0 • 

506 

WA  = 

PklM  AX  + 

BBC 

( 1)*AL0G(DELTAR/RR1) 

WC- 

0. 

DO 

66  N  =  2 

.M 

XN* 

N 

WB  = 

BBC ( N  ) 

* ALOG ( XN/( XN-1. 

)  ) 

WC- 

WC  +  WB 

66 

CON 

TINUE 

WD  = 

WA  +  WC 

W(  I 

XY)  =( 

.00 

708*H 

*DPR ) / WD 

IF 

(HIXY. 

NE. 

GIXY) 

GO  TO 

50 
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PLE  NO  1  TORTRAN  SOURCE  LIST  NINAL 

SOURCE  STATEMENT 

WRITE  (6,5004)  W(IXY) 

WRITE  (6,5005)  ( BBC ( 1 J ) , I J= I , M  ) 

50  RETURN 
END 


jahih  i2ij  aowuoe  mah^oi 


I  OKI  3J«1 

1 TAT2  33SU02 


lH,!*L!»UiU«fl|  t200e*d>  311** 
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15C052  EXAMPLE  NO  1 


IBLDR  —  JOB 


LEAL 


CT  PROGRAM  IS  BEING  ENTERED  INTO  STORAGE. 


63.000 
0.645 
0.230 
0.087 
X  .  700 


131. 5C0 
52.280 
0.031 
400.000 
1.200 


80.750 
5  1.450 
1.090 
0.266 
612.000 


71.900 

140.000 

530.000 

25.000 


0.260 
1.030 
95.000 
6.  120 


0.560 

0.500 

0.913 

14.392 


ORIGINAL  SATURATION 

0.200  0.200  0.200  0.200  0.200  0.200  0.200  0.200  0.200  0.540 
TOTAL  TIME  OF  INJECTION 

720.0  710.0  690.0  660.0  620.0  560.0  485.0  400.0  285.0  145.0 


INITIAL  BACK  PRODUCTION  RATE* 


384.31  BPD 


RESERVOIR  AFTER  4.0  DAYS 

TEMPERATURE  PROFILE 
526.6  524.5  524.8  524.4  524.3 
SOLUTION  ERROR 

-0.000  0.000  0.000  0.000  0.000 

RATE  367.3 

ELEMENT  BY  ELEMENT  RESISTANCE 

0.000  0.000  0.000  0.000  0.000 


524.1  523.7  522.9  519.3  478.6 

0.000  0.000-0.000  0.000-0.000 

0.000  0.000  0.000  0.000  0.004 


RESERVOIR  AFTER  8.0  DAYS 

TEMPERATURE  PROFILE 

462.9  424.3  427.2  421.7  417.6  410.1  399.4  384.0  353.4  275.1 
SOLUTION  ERROR 

-0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000 

RATE  359.9 

ELEMENT  BY  ELEMENT  RESISTANCE 

0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.001  0.007 
INCREMENTAL  AND  TOTAL  PRODUCTION  1454.3  4494.7 

SATURATION  DISTRIBUTION 

0.200  0.200  0.200  0.200  0.200  0.200  0.200  0.490  0.490  0.540 


RESERVOIR  AFTER  12.0  DAYS 

TEMPERATURE  PROFILE 

381.1  366.5  370.6  363.9  361.2  354.8  346.2  333.5  308.6  249.4 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000  0.000  0.000-0.000  0.000 

RATE  355.2 

ELEMENT  BY  ELEMENT  RESISTANCE 

0.000  0.000  0.000  0.000  0.000  0.000  0.001  0.003  0.003  0.008 
INCREMENTAL  AND  TOTAL  PRODUCTION  1430.1  5924.9 

SATURATION  DISTRIBUTION 

0.200  0.200  0.200  0.200  0.200  0.200  0.490  0.490  0.490  0.690 


95.0 


95.0 


95.0 


OdS.O 

OSI  ,d 


.  3t  !>■  .  2 


0  OOS.O  OOS.O  GOS.O  OQS.O  GOS.<  CM mO  os  5  ooslo 


ddS.O  000. 00* 

'00. SI*  00$,  f 


yfoi  T^SM  Mi  TO  3MIT  JA TOT 
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RESERVOIR  AFTER  16.0  DAYS 

TEMPERATURE  PROFILE 

354.1  350.2  350.4  344.4  341.4  334.1  324.5 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  346.6 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.000  0.000  0.000  0.000  0.000  0.001  0.003 
INCREMENTAL  AND  TOTAL  PRODUCTION  1403.6 
SATURATION  DISTRIBUTION 
0.200  0.200  0.200  0.200  0.200  0.490  0.490 


310.5  283.9  224.0 
0.000  0.000  0.000 

0.003  0.004  0.032 
7328.5 

0.490  0.690  0.690 


95.0 


RESERVOIR  AFTER  20.0  DAYS 

TEMPERATURE  PROFILE 

338.9  331.6  331.7  326.2  323.0  315.6  306.0 
SOLUTION  ERROR 

-0. 000-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  339.5 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.000  0.000  0.000  0.000  0.001  0.003  0.003 
INCREMENTAL  AND  TOTAL  PRODUCTION  1372.3 
SATURATION  DISTRIBUTION 
0.200  0.200  0.490  0.200  0.200  0.490  0.490 


292.0  265.7  209.6 
0.000  0.000  0.000 

0.003  0.017  0.034 
8700.8 

0.490  0.690  0.690 


95.0 


RESERVOIR  AFTER  24.0  DAYS 

TEMPERATURE  PROFILE 

323.1  318.5  318.3  312.8  309.5  302.0  292.3 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  335.6 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.000  0.000  0.000  0.000  0.001  0.003  0.003 
INCREMENTAL  AND  TOTAL  PRODUCTION  1350.3 
SATURATION  DISTRIBUTION 
0.490  0.490  0.490  0.490  0.490  0.490  0.490 


278.2  252.1  198.2 
0.000  0.000  0.000 

0.003  0.018  0.037 
10051.1 

0.690  0.690  0.690 


95.0 


RESERVOIR  AFTER  28.0  DAYS 

TEMPERATURE  PROFILE 

311.6  307.0  306.5  301.3  297.9  290.4  280.7 
SOLUTION  ERROR 

-0.000-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  309.8 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.003  0.003  0.003  0.003  0.003  0.003  0.003 
INCREMENTAL  AND  TOTAL  PRODUCTION  1290.8 
SATURATION  DISTRIBUTION 
0.490  0.490  0.490  0.490  0.490  0.490  0.690 


266.7  241.2  189.9 
0.000  0.000-0.000 

0.015  0.020  0.038 
11341.8 

0.690  0.690  0.690 


95.0 
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RESERVOIR  AFTER  60.0  DAYS 

TEMPERATURE  PROFILE 


257.7  254.8  253.8 
SOLUTION  ERROR 

249.6 

246.  1 

239.3 

230.6 

218.2 

196.7 

157.6 

-0.000  0.000  0.000 
RATE  223.3 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.015  0.015  0.015 

0.015 

0.016 

0.017 

0.018 

0.021 

0.027 

0.105 

INCREMENTAL  AND  TOTAL  PRODUCTION 

899.9 

19725.7 

SATURATION  DISTRIBUTION 

0.690  0-690  0.690 

0.690 

0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

RESERVOIR  AFTER 
TEMPERATURE  PROFILE 

64.0 

DAYS 

253.5  250.8  249.6 
SOLUTION  ERROR 

245.6 

242.  1 

235.4 

226.8 

214.6 

193.4 

155.4 

0.000  0.000  0.000 
RATE  220.1 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.015  0.015  0.015 

0.016 

0.017 

0.017 

0.019 

0.022 

0.028 

0.108 

INCREMENTAL  AND  TOTAL  PRODUCTION 

886.8 

20612.5 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

RESERVOIR  AFTER  68.0  DAYS 

TEMPERATURE  PROFILE 

249.5  247.0  245.8  241.8  238.4 

231.8 

223.3 

211.2 

190.5 

153.3 

SOLUTION  ERROR 

-0.000-0.000  0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

RATE  217.0 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.015  0.016  0.016  0.016  0.017 

0.018 

0.019 

0.022 

0.029 

0.110 

INCREMENTAL  AND  TOTAL  PRODUCTION 

874.2 

21486.7 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

RESERVOIR  AFTER  72.0  DAYS 

TEMPERATURE  PROFILE 

245.8  243.4  242.3  238.4  235.0 

228.5 

220.1 

208.2 

187.8 

151.5 

SOLUTION  ERROR 

-0.000-0.000  0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

RATE  214.1 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.016  0.016  0.016  0.017  0.017 

0.018 

0.020 

0.023 

0.029 

0.112 

INCREMENTAL  AND  TOTAL  PRODUCTION 

862.1 

22348.9 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

95.0 


95.0 


95.0 


95.0 


000.0 

80i  .0 

000.0 

01  0 
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RESERVOIR  AFTER 
TEMPERATURE  PROFILE 

80.0 

DAYS 

239.2  237.1  235.9 
SOLUTION  ERROR 

232.  1 

228. 8 

222.5 

214.3 

202.  7 

182.9 

146. 1 

0.000  0.000  0.000 
RATE  208.4 

0.000 

0.000 

0.000- 

O.OCO 

0.000 

0.000 

0.000 

ELEMENT  8Y  ELEMENT 

RESISTANCE 

0.017  0.017  0.017 

0.018 

0.018 

0.019 

0.021 

0.024 

0.030 

0.115 

INCREMENTAL  AND  TOTAL  PRODUCTION 

839.2 

240 

38.6 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

RESERVOIR  AFTER  84.0  DAYS 

TEMPERATURE  PROFILE 


236.2  234.2  233.0 

229.3  226.0 

219.7 

211.7 

200.2 

180.7 

146.6 

SOLUTION  ERROR 
-0.000-0.000  0.000 

0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

RATE  205.7 

ELEMENT  BY  ELEMENT 
0.017  0.017  0.017 

RESISTANCE 
0.018  0.019 

0.020 

0.021 

0.024 

0.031 

0.117 

INCREMENTAL  AND  TOTAL  PRODUCTION 

828.2 

24866.8 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

RESERVOIR  AFTER  88.0  DAYS 

TEMPERATURE  PROFILE 


233.4  231.5  230.3 

226.6  223.3 

217.2 

209.2 

197.9 

178.7 

145.3 

SOLUTION  ERROR 
0.000-0.000  0.000 

0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

RATE  203.1 

ELEMENT  BY  ELEMENT 
0.017  0.018  0.018 

RESISTANCE 
0.018  0.019 

0.020 

0.022 

0.025 

0.031 

0.119 

INCREMENTAL  AND  TOTAL  PRODUCTION 

817.6 

25684.4 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

RESERVOIR  AFTER  92.0  DAYS 

TEMPERATURE  PROFILE 


230.7  228.9  227.7 

224.1  220.9 

214.8 

206.9 

195.7 

176.8 

144.0 

SOLUTION  ERROR 
-0.000-0.000  0.000 

0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

RATE  200.5 

ELEMENT  BY  ELEMENT 
0.018  0.018  0.018 

RESISTANCE 
0.019  0.019 

0.020 

0.022 

0.025 

0.032 

0.120 

INCREMENTAL  AND  TOTAL  PRODUCTION 

807.1 

26491.5 

SATURATION  DISTRIBUTION 

0.690  0.690  0.690  0.690  0.690 

0.690 

0.690 

0.690 

0.690 

0.820 

95.0 


95.0 


95.0 


95.0 
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RESERVOIR  AFTER  572.0  DAYS 

TEMPERATURE  PROFILE 
152.1  151.8  151.0  149.5  147.6  144.5 
SOLUTION  ERROR 

-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  106.9 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.068  0.068  0.070  0.072  0.075  0.079 
INCREMENTAL  AND  TOTAL  PRODUCTION 
SATURATION  DISTRIBUTION 
0.820  0.820  0.820  0.820  0.820  0.820 


RESERVOIR  AFTER  576.0  DAYS 

TEMPERATURE  PROFILE 
151.9  151.6  150.8  149.3  147.4  144.4 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000 

RATE  106.7 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.068  0.069  0.070  0.072  0.075  0.079 
INCREMENTAL  AND  TOTAL  PRODUCTION 
SATURATION  DISTRIBUTION 
0.820  0.820  0.820  0.820  0.820  0.820 


RESERVOIR  AFTER  580.0  DAYS 

TEMPERATURE  PROFILE 
151.7  151.4  150.6  149.1  147.2  144.2 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000 

RATE  106.5 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.068  0.069  0.070  0.072  0.075  0.080 
INCREMENTAL  AND  ..TOTAL  PRODUCTION 
SATURATION  DISTRIBUTION 
0.820  0.820  0.820  0.820  0.820  0.820 


RESERVOIR  AFTER  584.0  DAYS 

TEMPERATURE  PROFILE 
151.5  151.2  150.4  148.9  147.0  144.0 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000 

RATE  106.3 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.068  0.069  0.070  0.072  0.075  0.080 
INCREMENTAL  AND  TOTAL  PRODUCTION 
SATURATION  DISTRIBUTION 
0.820  0.820  0.820  0.820  0.820  0.820 
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134.5 

125.3 

111.6 

95.0 

0.000 

0.000 

0.000 

0.000 

0.086 

428.0 

0.099  0.122 
89665.2 

0.168 

0.820 

0.820 

0.820 

0.820 

140.3 

134.4 

125.2 

111.6 

95.0 

O.OCO 

0.000 

OiOOO 

0.000 

0.087 

427.2 

0.099  0.123 
90092.4 

0.168 

0.820 

0.820 

0.820 

0.820 

140.1 

134.2 

125.1 

111.5 

95.0 

0.000 
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0.000 

0.000 
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426.4 
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0.820 
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150052  EXAMPLE  NO  2 

IBL  OR 

—  JOB 

LEAL 

ECT  PROGRAM  IS  BEING 

ENTERED 

INTO  STORAGE. 

42.790 

1 18. 8C0 

61.730  52.250 

0.2  30 

0.2  60 

0.545 

60. 1 30 

50.320  140.000 

1.013 

0.500 

0.230 

0.042 

1.500  445.000 

74.000 

0.7  65 

0.235 

1600.000 

0.362  10.000 

0.945 

54.945 

1.700 

1.200 

3000.000 

ORIGINAL  SATURATION 

0.200  0.200 

0.200  0 

.200  0.200  0.200  0.200 

0.200  0. 

370  0.370 

TOTAL  TIME  OF  INJECTION 

720.0  710.0  690.0  660.0  620.0  560.0  485.0  400.0  285.0  145.0 


INITIAL  BACK  PRODUCTION  R ATE  = 


103.64  BPD 


RESERVOIR  AFTER 
TEMPERATURE  PROFILE 


4.0  DAYS 


436.7  432.6  433.2 
SOLUTION  ERROR 

432.4 

432.  2 

431.6 

430.7 

429.0 

422.7 

375.4 

0.000  0.000  0.000 
RATE  103.3 

0.000 

0.000 

o.ooo- 

O.OCO 

0. 000- 

■0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

RESERVUIR  AFTER 
TEMPERATURE  PROFILE 

8.0 

DAYS 

373.9  347.5  349.5 

344.4 

340.7 

3  3  3.7 

323.7 

309.  1 

279.2 

207.2 

SOLUTIGN  ERROR 
-0.000-0.000  0.000 
RATE  100.7 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

INCREMENTAL  AND  TOTAL  PRODUCTION 

408.2 

1236. 7 

0.200  0.200 

0.200 

0. 

200 

0.200 

0.200 

0.435  0.200 

0.370 

0. 

630 

RESERVOIR  AFTER 
TEMPERATURE  PROFILE 

12.0 

DAYS 

307.1  301.0  303.0 

297.6 

295.2 

289.5 

281.8 

270.0 

246.4 

192.9 

SOLUTION  ERROR 
-0.000-0.000  0.000 
RATE  99.4 

0.000 

0.000 

0.000 

0.000 

0.  000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.008 

INCREMENTAL  AND  TOTAL  PRODUCTION 

400.2 

1636. 9 

0.200 

0.435 


0.200 

0.2C0 


0.200 

0.630 


0.200 

0.630 


0.200 


74.0 


74.0 


74.0 


0.435 
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RESERVOIR  AETER  16.0  DAYS 

TEMPERATURE  PROFILE 
290.0  287.5  286.5  282.0  279.0  272.5 
SOLUTION  ERROR 

-0.000-0.000  0.000  0.000  0.000  0.000 

RATE  98.1 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.000  0.000  0.000  0.000  0.000  0.000 
INCREMENTAL  AND  TOTAL  PRODUCTION 

0.200  0.200  0.200  0 

0.435  0.200  0.630  0 


263.8  250.7  225.3 

0.000  0.000  0.000 


0.000  0.000  0.004 
394.9  2031.8 

200  0.200 

630 


170.9 

0.000 


0.009 

0.435 


RESERVOIR  AFTER  20.0  DAYS 

TEMPERATURE  PROFILE 

276.2  271.4  271.0  266.7  263.6  257.1  248.4 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  96.8 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.000  0.000  0.000  0.000  0.000  0.000  0.000 
INCREMENTAL  AND  TOTAL  PRODUCTION  389.7 

0.200  0.2C0  0.200  0.435 

0.435  0.630  0.630  0.630 


235.3  210.4 

0.000  0.000 


0.000  0.005 
2421.5 
0.435 


160.5 

0.000 


0.009 

0.435 


RESERVOIR  AFTER 
TEMPERATURE  PROFILE 


24.0  DAYS 


262.9  260.6  259.7 
SOLUTION  ERROR 

255.5 

252.3 

245.7 

236.9 

223.7 

199.2 

151.4 

0.000-0.000  0.000 
RATE  95.5 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.000  0.000  0.000 

0.000 

0.000 

0.000 

0.001 

0.004 

0.005 

0.010 

INCREMENTAL  AND  TOTAL  PRODUCTION 

384.6 

2806.2 

0.435  0.435 

0.435 

0. 

435 

0.435 

0.435 

0.630  0.630 

0.630 

0. 

630 

RESERVUIR  AFTER 

28.0 

DAYS 

TEMPERATURE  PROFILE 

253.8  251.1  250.0 
SOLUTION  ERROR 

246.0 

242.8 

236.2 

227.4 

214.2 

190.1 

144.7 

0.000-0.000  0.000 
RATE  94.2 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RES  I  STANCE 

0.000  0.000  0.000 

0.000 

0.000 

0.001 

0.004 

0.004 

0.006 

0.010 

INCREMENTAL 
0.435 
0.6  30 


AND  TOTAL 
0.435 
0.630 


PRODUCTION 

0.435 

0.630 


379.5 

0.435 

0.630 


3185.6 

0.435 


0.435 


74.0 


74.0 


74.0 


74.0 
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RESERVOIR  AFTER  120.0  DAYS 

TEMPERATURE  PROFILE 


173.3  172.7  171.4  16B.9 

165.  9 

160.7 

153.9 

143.8 

127.5 

102.9 

SOLUTION  ERROR 
-0.000-0.000  0.000  0.000 
RATE  67.2 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT  RESISTANCE 

0.005  0.005  0.006  0.006 

0.006 

0.006 

0.007 

0.008 

0.009 

0.013 

INCREMENTAL  AND  TOTAL  PRODUCTION 

271.0 

10535.5 

0.630  0.630 

0.630 

0. 

630 

0.630 

0.630 

0.630  0.630 

0.630 

0. 

630 

RESERVOIR  AFTER  124.0 

TEMPERATURE  PROFILE 

DAYS 

171. R  171.2  169.9  167.5 
SOLUTION  ERROR 

164.  5 

159.4 

152.6 

142.6 

126.5 

102.3 

0.000-0.000  0.000  0.000 
RATE  66.1 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT  RESISTANCE 

0.006  0.006  0.006  0.006 

0.006 

0.006 

0.007 

0.008 

0.009 

0.013 

INCREMENTAL  AND  TOTAL  PRODUCTION 

266.7 

10802.2 

0.630  0.630 

0.630 

0. 

630 

0.630 

0.630 

0.630  0.630 

0.630 

0. 

630 

RESERVOIR  AFTER  128.0 

TEMPERATURE  PROFILE 

DAYS 

* 

170.4  169.8  168.6  166.1 
SOLUTION  ERROR 

163.  1 

158.  i 

151.4 

141.4 

125.5 

101.7 

0.000-0.000  0.000  0.000 
RATE  65.0 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT  RESISTANCE 

0.006  0.006  0.006  0.006 

0.006 

0.006 

0.007 

0.008 

0.009 

0.013 

INCREMENTAL  AND  TOTAL  PRODUCTION  2 62.3 

0.630  0.630  0.630  0.630 

0.630  0.630  0.630  0.630 


11064.5 

0.630 


0.630 


74.0 


74.0 


74.0 


RESERVOIR  AFTER  132.0  DAYS 

TEMPERATURE  PROFILE 
169.0  168.5  167.2  164.8  161.' 
SOLUTION  ERROR 

0.000-0.000  0.000  0.000  0.00 

RATE  64.0 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.006  0.006  0.006  0.006  0.00 
INCREMENTAL  AND  TOTAL  PRODUCTION 
0.630  0.630  0.630 

0.630  0.630  0.630 


156.9 

150.2 

140.3  124.6 

101.2 

0.000 

0.000 

0.000  0.000 

0.000 

0.006 

0.007 

0.008  0.010 

0.013 

ON 

258.0 

11322.5 

0. 

630 

0.630 

0.630 

74.0 


0.630 
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RESERVOIR  AFTER 
TEMPERATURE  PROFILE 

360.0 

DAYS 

132.9  132.7  131.8 

130.2 

128.0 

124.6 

119.9 

113.0 

102.9 

89.  1 

SOLUTION  ERROR 
-0.000-0.000  0.000 
RATE  42.6 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

ELEMENT  BY  ELEMENT 

RESISTANCE 

0.102  0.104  0.108 

0.115 

0. 127 

0.  146 

0.179 

0.247 

0.408 

0.818 

INCREMENTAL  AND  TOTAL  PRODUCTION 

170.5 

22470.8 

0.830  0.830 

0.830 

0. 

830 

0.830 

0.830 

0.830 


0.830 


0.830 


0.830 


74.0 


RESERVOIR  AFTER  364.0  DAYS 

TEMPERATURE  PROFILE 
132.6  132.4  131.5  129.9  127.7 
SOLUTION  ERROR 

-0.000-0.000  0.000  0.000  0.000 

RATE  42.5 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.103  0.105  0.110  0.117  0.128 
INCREMENTAL  AND  TOTAL  PRODUCTI 


0.830 

0.830 


0.830 

0.830 


0.830 

0.830 


124.3 

119.6 

112.8  102.7 

89.0 

0.000 

0.000 

0.000  0.000 

0.000 

0.147 

0.181 

0.249  0.411 

0.820 

ON 

170.1 

22640.9 

0. 

830 

0.830 

0.830 

74.0 


0.830 


RESERVOIR  AFTER  368.0  DAYS 

TEMPERATURE  PROFILE 
132.3  132.1  131.2  129.6  127.4  124.0 
SOLUTION  ERROR 

-0.000  0.000  0.000  0.000  0.000  0.000 

RATE  42.4 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.104  0.106  0.111  0.118  0.129  0.149 
INCREMENTAL  AND  TOTAL  PRODUCTION 

0.830  0.830  0.830  0 

0.830  0.830  0.830  0.830 


119.3 

112.6 

102.5 

88.9 

0.000 

0.000 

0.000 

0.000 

0.182 

0.251 

0.413 

0.822 

169.7 

228 

10.6 

830 

0.830 

0.830 

74.0 


RESERVOIR  AFTER  372.0  DAYS 

TEMPERATURE  PROFILE 
131.9  131.7  130.9  129.3  127.1 
SOLUTION  ERROR 

-0.000-0.000  0.000  0.000  0.000 

RATE  42.3 

ELEMENT  BY  ELEMENT  RESISTANCE 
0.106  0.107  0.112  0.119  0.131 
INCREMENTAL  AND  TOTAL  PRODUCTI 
0.830  0.830  0.830 

0.830  0.830  0.830 


123.7 

119.1 

112.3  102.4 

88.8 

0.000 

0.000 

0.000  0.000 

0.000 

0.  150 

0.184 

0.253  0.416 

0.824 

ON 

169.3 

22979.9 

0. 

830 

0.830 

0.830 

74.0 


0.830 
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